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NOMENCLATURE 

The following is the list of some important symbols 
and notations used in this work. 

u(z,s) = out-of -plane displacement i.e, longitudinal 

displacement in z- direction, 
v( 2 ,s) = in-plane tangential displacement, 

w(z,s) = in-plane normal displacement. 

u_ = uniform longitudinal extention of the cross-section 

Q, 

* displacement of a point B(Fig,2) along 

X y 

X and y axes 

“ rotations at the point B about x and y axes 

X y 

= angle of twist of the cross-section 

^ = qL and is a measure of warping 

^b 

x,y * coordinates of point A 

p,q = perpendicular distances from the point B to 

the tangent and normal at point A respectively . 

= warping function (sectorial coordinate) . 

t thickness of the profile 

? “ torsional function 

f shear flow 

G » shear ie^uIus 



(O 


differentiation wrt z-coordinate 


(-) 

T 

U 


0 

z 

■^zs 

(.) 

p 

E 

A 


XX 


"(«)W 


rr 


PP 


qq 




rq 


pq 


differentiation wrt s -coordinate 

kinetic energy of the beam 
strain energy of the beam 

normal stress 

shear stress 

differentiation wrt time 

mass density of the beam 
modulus of elasticity 

area of the cross-section 

second moment of area about y-axis 

second moment of area about x-axis 

sectorial moment of inertia 

torsional constant 

warping rotational moment of inertia 
shear area about x-axis 
shear area about y-axis 
shear area moi^nt about x-axis 
shear area moment ahcmt y-axis 
mixed shear area 
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polar moment of inertia of the cross section 
about the center of shear 


(e) 

I J 

(ne) 

[ ] 

T 

C ] 

IN j 

|F}(ne) ^ 

and q^, 

®y 


variation 

denotes the typical element 

row vector 
column vector 

at the nodes of typical element 

square matrix 

transpose of matrix 

interpolation functions 

mass matrix of the element 

stiffness matrix of the element 

displacement vector of the element 

boundary force vector of the element 

work done by external- forces 

line loads per unit length of the beam in the 
x,y and z directions respectively 

applied torsional moment ber»iing moment 
(about x), bending moment (about y) atisd 
bimoments per unit length of beam , 
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and = 


concentrated load in the x,y and z directions 
respectively 


Mt,Mx,My = 

and 


concentrated torsional torsional moment, 
bending moment (about x), bending moment 
(about y) , and bimoments. 


6(z-z) 


U f u 


Dirac delta function for concentrated loads 
and moments . 

load vector of the element 

the displacements of pole B in the x and 
y-directions respectively 




coordinates of the pole B 

angle of tangent at A with the x-axis 


the perpendicular distance from the point 0 to 
the targent at point A 


^o 


the perpendicular distance from the point 0 to 
the normal at point A 


(Xq.yq) = 

[T] 

Wxx 

#(Sg) 


coordinates of the point A 

coordinates of the point of loading 

transformation matrix 
bending moB^nt 
bimoment 

warping function 
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^st. 


shear force 

statical moment about x-axis 
thickness of the profile 

warping twisting moment or flexural twist 
statical sectorial moment 
Saint Venant twisting moment 

Saint venant torsional constant of the section 
warping normal stress 
warping shear stress 
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SYNOPSIS 

In this work free vibration, static and buckling 
analyses of thin-walled beams are done. The effects of 
warping, rotary inertia and shear deformation on the frequencies 
and displacements are studied. Effect of warping on the 
critical buckling loads is studied, Monosymmetric and unsymmetric 
cross sectional beams are considered. Torsional analysis of a 
box section and tapered I-section beams is also done. Plane 
frames are analysed for static and free vibration cases. 

The governing differential equations are derived from 
the variational principles, and are non-dimensionalised. The 
finite element method is used to solve these equations , The 
field variables are approximated over the elements by polynomials 
The residues obtained are minimized by the Galerkin method. 

Gauss quadrature is used for obtaining the elemental matrices. 

The elemental matrices are assembled to get the global matrices 
and these are solved for frequencies, displacements and 
critical buckling loads, for various boundary conditions. 
Equations of the plane frames in the ^obal coordinates are 
obtained by using transformation matrices. 
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While considering shear effects, FEM solutions are 
obtained from six coupled second order differential equations 
as this enabled to satisfy the geometrical and natural 
boundary conditions .dhree coupled fourth order differential 
equations are used when the shear effects are neglected. 

The effect of warping, rotary inertia, and shear on 
the frequencies of vibration is studied. A little parametric 
study is done with the non -dimens ionalised parameters. In the 
case of tapered beams the effect of taper on the displacements 
is studied , 

In the case of plane frames the effect of warping and 
configuration of the beams on the displacements and stresses is 
studied. Frequencies of vibration of plane frames are obtained, 
when the effects of warping and rotary inertia are taken into 
account , 



CHAPl'Eii- I 


IlfrRODUCTION 

It is very difficult to define precisely a thin walled 
beam. One may say that it is an assembly of thin plates joined 
alon£, their edges. The plate thicknesses are small compared 
to the cross-sectional dimensions, and cross sectional 
dimensions inturn are small compared to the length of the 
beam. Thin walled beams have lot of applications in the field 
of engineering owing to their high strenth to low weight 
ratio. 

The major application of the thin-v/alled beams is in 
the aircraft industry, ship structures, cranes, bridge 
construction, storage racks, railway wagons and coaches, 
Thin-valled beams can be used for different purposes depending 
upon the design requirements. The beams can be used as closed 
section box girders to give more torsional rigidity and as 
open section plate girders to give low torsional rigidity . 

In thin-v/alled beams the shear stresses and strains 
are relatively much larger than those in solid beams. When a 
thin-v;alled beam is tv/isted, there is an out of plane 
distortion of the cross section of the beam in the dii*ectioa 
of the beam axis, called warping. This violates the Bernoulli* 
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hypothesis (plane sections remain plane over the entire 
cross section) . If the ..'arpin^ of the cross section is restrains 
it results in producin.:^ normal stresses in the axial direction 
and shear stresses in the cross section. These stresses are 
called the v/arpin^' normal stress (a^) and warping shear stress 
(t^) or in short warping stresses. These .■■az-ping stresses are 
as larjj^e as or even larger than the bending stresses and 
can not bo ignczred. Thus the thin-wailed beams have to be 
treated in a special manner and a separate theory has been 
developed to analyse these beams^. i .e Vlasov’ s theory for the ti 
walled beams. It is more general theory than the Bernoulli’s 
theory and is based on the law of sectorial areas and treats 
the plane sections case as a special case. Two new concepts, 
bimoment and v/arping tv/isting moment, were introduced by 
Vlasov’s theory to explain the warping phenomena and the stress« 
developed there in. 

In the case of thin-v;alled beams there can be a local 
buckling, if the inplane stresses reach their critical values. 
This changes the geometry of the cross section of the beam 
in contrast to the over all buckling where the geometry of tbe 
cross section is retained. If the thin-walled beam is sufficieni 
long then over all buctiing occurs before local buckili^ 
can occur. In some cases these two buoklings^ can interact and 
in tjiose cases the buckling load will be loss than the individu? 
buclfling loads. Thus the thin-walled beams have to be designed 
both overall Aiad Inoal hiioldi 1 no-a . 
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There ie 'm irripor bent phenomena oi' coupled flexure 
and torsion which occurs in the thin-walled beams of mono and 
unsymmetric cross section be.'ms. This phenomena can be explained 
v/ith respect to the she-'^r center. Sheer center is a point 
through v/hich the transverse sneer loads and supports should 
act, rather than througn the centroid of the cross section, 
if no torsion is to result. Hence if one shear center does not 
coincide with, the centroid of the cro 'S section (e.g. mono- 
symmetric and unsyiratetric cross sections) and if the transverse 
loads act through the centroid, then there exists a coupling 
betwee/i flexure and torsion. 

1.1 Pr.£VIOU3 'u-ORK : 

The coupled flexure and torsion, 'which is an 
interesting phenomena of thin-v;alled beams, has been studied 
by many workers and researchers in this field . 

Lot of work was done on the static analysis of thin- 
'walled beams. Number of books are available in the literature. 
Vlasov [i ] in his very comprehensive famous book 
"*Thin-.valled. Elastic Beams" developed his theory far static, 
stability and dynamics for thin-wallod straight and curved 
beams. Among the some other important books dealing with 
static analyses are ?^urray [2 ] , Goelsvih. [3 ], 

Megs bn [4 ] snd Zbirohovslii [5 ]. 
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In hi3 booK, [2 ] 'erivcd thi'' ecjuilibrium 

equations for bendinn. and torsion of bhin-v.'olled beasns by first 
order tnoory. The e';ur Liens vjer>^ derived by both equilibrium 
and energy inufchoda. Some e/oamples on torsional analysis of 
thin-v;alled beams ''.'ere solved. 

O^elsvil: [3 ] discussed the torsion and flexure of 
thin-v.'uiled bars of open and closed cross sections. The 
section properties and imj-ortanb features of regular sections 
V/'ere discussed . 

The books by Megson [4 ] and Zbirohowski [5 ] 
also discussed mariy aspects of thin-walled beams, under static 
loading . 

A finite element fonnulation for thin-walled beams was 
given by Gunnlaugsson [ 6 ] • this work the governing 
equilibrium equations v/ero drived by including the deformation 
due to shear. A transformation matrix w'as given in the vvork 
to change the local no al parameters to tlie global nodal 
parameters in order to vochieve the compatibility betv/een the 
tv/o adjacent elements, whose centroidal axes and. axes of twist 
do not coincide. 

A torsional stiffness matrix was derived by Chen 8urad 
Hu [ 7 ] » a thin-ualled beam b ;sed on the warping 
diSiAacement assumption. Also the general stiffness matrix 
for a beam under simultaneous torsion, bending and tension 
was derived for the coupled torsion-bending analysis . 
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L transformation matrix to chan^^e tl.C’ local coordimtes of an 
element to the global system v/as given in their work,. The 
correctness of this was checked for a thin walled beam of 
rectangSJlar box cross section under uniform distributed torque, 

V/ekezer [8 ] applied the finite element method to' 
the theory of thin-walled bars of variable cross sections in 
torsional analysis 1 The solution of the problem was based on 
the linear membrane shell theory with the application of 
Vlasovas assumptions. He divided the bar into elements along 
its longitudinal axis and then approximated the shell mid 
surface of the element by arbitrary triangular sub -elements . 

There has been little work in the dynamic analysis 
of thin-walled beams under coupled flexure and torsion, 

Vlasov [ 1 ] , did work on coupled flexural torsional vibration: 

of thin walled beams including the effects of warping and 
rotary inertia. The applicability of this coupled flexure and 
torsion to many practical problems was considered by many 
researchers using analytical as ^.'ell as numv rical methods. 

The tripled coupled flexural and torsional vibrations of 
thin-waHiEKi beams were derived by Gere and Lin [ 9 ] based on 
the development given by Timoshenko [10 ] for dcwble coupling. 
The effects of rotary inertia and shear deformation were 
neglected. The effect of warping was taken into account. 
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Bishop and Price [15] solved bhe enuntion of coupled 
bending and twisting of a Timoshenko beam. The effects of 
shear deflection and rotcry inertia in bending for an 
unsymmetric beam were included. 'The mode shapes and the 
orthogonality conditions were obtained. 

The effects of shear deformo. tion, rotary inertia and 
v/arping were considered by Tso [16] , in the coupled flexural 
torsional vibrations of thin-walled be.-.iras . The governing equations 
of motion v/ere derived by variational principles. The coupled 
frequencies of vibrations of a raonosymmetric cross section 
beam were obtained by analytical methods. Tso discussed the 
effect of shear and rotary inertia on the frequencies , 

The coupled flexural torsional buckling of thin walled 
columns has been considered by many workers. The books dealing 
this topic are Timoshenko and Gere {17 "] , Vlasov [i ] 

and Murray [2 ] . These obtained the governing equations for 

the coupled flexural torsional buckling "by the equilibrium * 
method. The equations for monosymmetric and unsymmetric 
section beams were solved by exact methods. Critical load 
equations for monosymmetric and unsymmetric section beams, 
with any end conditions v;ere given, '.ts functions of uncoupled 
critical buckling loads in flexure and torsi on ♦ 

The analysis of a plane frame v/as considered by many • ■ 

workers in the past, ffertin [18 ]> gave the transformation 
matrices for Euler beam elements, in order to change the local 


1 "nn -irVv/rv rrl 
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solved by direct stiffness and flexibility methods, 

Nosseir and Dickinson [19] applied the finite element 
method to the free vibration analysis of a model of a car body, 
^hey approximated the chassis, under body and full body of a 
car with ordinary beam elements and panels , The out of plane 
•vibrations were considered. Because oC the symmetry of the 
frame, only one half of the model v;ao analysed with appropriate 
boundary conditions along the axis of symmetry. 

Finite element techniques were applied by Ali, Hedges, 
Mills, Horville and Gurdoganj| 20] 'to the analysis of an 
automobile structure. The .static, stress and dynamic analyses 
v/ere done on an automobile chassis frsme. Ordinary beam elements 
were used in the approximation of the frame . 

1 .2 Pl-ldSSriT WORK : 

In the present work the Galerkin finite element method 
is used to solve the more familiar differential equation of 
the thin-walled beams and frames. First the equations of 
motion of a thin-walled beam of unsymraetric cross section 
with variable area were derived by variational principles, 
including the effects of v/arping, rotary inertia and shear for 
completeness sake. The effect of these quantities on the 
coupled frequencies of vibration Is studied, for various end 



conditions. Non-dimensionalisation is uono for the governing 
equations and a limited parametric study is -’one v.'ith respect 
to these non-dimensionalised parameters. Results for both 
monosymmetric and unsyrametric cross section booms are obtained 
in this work. 

Static equilibrium equations Cor a general unsymmetric 
section of a thin walled beam are derived for completeness 
sake, including the effects of v/nrpin,. and shear. The coupled 
equations are solved for various end conditions, to get the 
coupled displacements and rotations of thin-v/alled beams of 
monosymmetric as well as unsyrametric sections. Toi^sional 
analysis of a cantilever of a rectangular box section and a 
tapered cantilever of I -section is done. In the case of tapered 
cantilever, the equations of equilibrium derived for a variable 
area are used. 

Coupled flexural torsional buckling equation for an 
unsymmetric section with variable area are derived by variational 
principles , The coupled equations aro solved for an 
unsymmetric section column with various end conditions. 

A transformation matrix is derived for transforming 
displacements in local coordinates to the displacements in 
global coordinates for using in the plane frame analyses. ,A 
plane frame is analysed for static and free vibration cases. 

Thus the principal objective of this work is to study 
the effects of rotary inertia and shear deformation in thin- 

walled beams and frames in static, free vibration and 
buckliriK cases. 


Cffi^FrER- II 


BASIC AND FEM EQUATIONS 


In this chapter governing differential equations 
of thin v/alled heams for the following cases are derived 
from the variational principle. 

Free vibrations 
Static analysis 
Buckling analysis 
Plane frame 

These equations are for variable cross-sections and are 
non dimensionalised , The equations account for warping, 
rotatory inertia and shear deformation. Their FEM equations 
are obtained using Galerkin method. 

2.1 FREE VIBMTIONS OF THIN VALLJSD BEAMS ; 

Equations of motion have been drived for a thin-walled 
beam as shown in Fig. 1. The origin of the axes is at the 
center of gravity (C.G) of the cross-sections and xy-axes 
are its principal axes. The z- axis is in the direction of the 


beam axis 
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The following are the assumptions made in deriving 
the Equations [ 2,6,16,21 ] 


1 . The deformations are small compared with the 
structural dimensions . 

2. The beam is straight before loading and there is no 
local buckling . 

3. The thickness of the cross section of the beam is 
small compared to the' cross sectional dimensions. 

4. wlien a thin-walled beam is bent and b^isted its 
cross section remains undeformed in its own plane, 
but rotates about the flexural axis (locus of the 
centres of shear) as a rigid body, 

2.1 .1 The in plane and out-of- plane displacements of a 
point on the cross-section- : 


Based on the assumptions made above, the displacement 
field for a point A on the cross-section (Fig. 2) is, 

Eqns. (A-34),[ 2,6] 


u(z,s) = Ug^- Oy x+ Q^y - 0Q 
v(z,s) = u^ X + Uy y + pQ^ 
w(z,s) = -u^ y + Uy X -f qO^ 

Here 


.... ( 1 ) 
.... ( 2 ) 
.... (3) 


u(z,s) = out- of- plane displacement i.e, longitudinal 
displacement in z -direction, 
v(z,s) = in.p>lane tangential disp3.acement , 
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W(2,S) 




Q ,Q 
x’ y 


05 


x,y 

p»q 


CO 

CO 

(0 

t 

If 

or 1' 

T 

G 


in-plane normal displacement, 
uniform longitudinal extention of the 
cross-section 

displacement of a point B (Fig, 2) along 
X and y axes , 

rotations at the point B about x and y axes, 


angle of twist of the cross-section, 

0^ and is a measure of warping ■ 

^b 

coordinates of point A 

perpendicular distances from the point B 
to the tangent and normal at point A 
respectively , 

warping function (sectorial coordinate) and is, 
Eqns. (A-31 ) , [2 ]. 

; (p- :^)ds for closed profiles 

° .... (4) 

s 

^ p ds for open profiles 
thickness of the profile 

torsional function and is, Eqns, (A-29)» [2] 


_T 

G e ' 

Shear flow 
shear modulo s 


• • # • 


(5) 
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In the above a prime ( ') denotes differentiation 
v^rt. z-coordinate and over dash (-) indicates differentiation 
wrt . s-Goordinate . 

2,1.2 The Kinetic and strain energies of the beam : 

In obtaining the expresions for kinetic and strain 
energies the following assumptions are made 

(i) The material of the beam is isotropic, homogeneous and 
linearly e lastic . 

(ii) The longitudinal stress a and the in-plane shear stress 

T__ are the only significant stresses in the beam. Other 
stresses and are negligible. 

These expressions and the corresponding fourth order 
differential equations were derived by ] ^or uniform 
cross-sections . Here second order differential' equations are 
derived for variable cross-sections. 

The kinetic energy of the beam can be ■written as 

^ .... ( 6 ) 

o A 

where p is the mass density of the beam. 

V/ith the above assumptions in the mind, the strain 
energy of the beam can be written as 

u = ^ ; ; (e(||)1g(||,. |g)^ dA dz 

O A 


• . . « 


( 7 ) 
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t'lTiere S is the modulus of elasticity and G is the modulus 
of rigidity . 

Substituting the expressions of u,v and w from the Eqns . (1) 
to (3) in the Eqns. (6) and (7), one gets 

r-U= |{P [/u^ dA+0^ /x^dA-2u^ 4^ / x dA + 0^ / y%A 

+ 0^ ^ OJ^dA - 20^ 0 yodA + 2u^ 4 / ydA 

A ^ A ^ ^ A 

-24„ ©„ fxy dA-2u^ 0 / oidA + 20 ^ / xoidA+u^ /(x^+y^)dA 

xy^ ^A y li ^ A 

+ Uy /(x^+y^)dA + 0^ ;(p^+q^)dA+24 ©2 ;(x-Xg)dA 

ri A A 

-2^ i [ (y-yB)<iA]] 

-E [u ' 2 ; dA+42 /x^dA-au s' ; xdA+dJ ; y^iA+i ^ 

^A ^A ^^A ^ A ^ 

-2©' 0' JywdA +2u^0' / ydA -2©'©' / xydA- 

-2u^ 0^ J wdA + 20^ 0^ S x(ixiA ] 

* A y A 

-G[ (u -0^)^ / x^<iA+(u'+@ / y^dA+(©'-0)^ / p%A 

+ 2(uJ^-0y) (Uy+0^) / X y dA + 2(u^+©^)(©2- 0 ) y p dA 

A 2 ^ " 

+ 2(0'-0)(u -9 ) / X p dA+ 0^ / ^ dA +2(u -0 )0 ; | x dA 

^ ^ y A A xy^-c 

+2(u^+0^)0 / ^ y dA+2(02,-0)0 (w+ ^ )dA 3 j’dz 

A 

...» ( 8 ) 

Dot (,) represents- differentiation wrt time 



As the origin of the principal axes, xy, is at the centroid 
of the cross-section, of the "beam, one has 


First 

moment 

of 

area about 

y-axis = 1 ^= 

S X 

A 

dA = 

0 

First 

moment 

of 

area about 

x-axis = I = 

/y 

A 

dA = 

0 

Product moment 

of inertia 


/ xy dA 

= 0 


A 

.... (9) 

If point B (Fig, ’>2) is centre of shear, then the 
following quantities are zero [3 ] , see Appendix II. 

Sectorial product area about x-axis = I,_ = / wy dA= 0 

Sectorial product area about y-axis = 1^^^ = / wx dA = 0 

( 10 ) 

Since in the calculation of the v/arping function (sectional 
coordinate) the starting point V (Fig, 2) is arbitrary, 
it can be chosen such that 

First moment of sectorial area about center of shear 

= I,= / codA = 0 .... (11) 

^ A 

By Leibnitz rule ^22 ] and using Eqns, (9) and (11), 
the following three integrals of Eqn. ( 8 ) become zero 



/ I X 

A ^ 

^|y 

A . 

/| 5 

A ^ 


dA = | 

^ dA - 


d 

3? 

d.A = | 

^ S dA = 

A ds 

1 

t 

d 

ds 

dA = | 

f dw , 

/ 35 

t 

. d 
3s 


i\ 

J y dA = 0 
A 

■f w dA = 0 
A 


... ( 12 ) 


Various quantities of Eqn. (8) are defined as follow 

Area of the cross-section = A = / dA 

A 

p 

Second moment of Area about y-axis= I„„= / x dA 

yy ^ 

Second moment of Area about x-axis = I = / dA 

XX ^ 

Sectorial moment of inertia = ^ u^dA 

2 

Torsional constant = Ij = / ~o dA 

A t^ 

Warping rotational moment of inertia= I / p^dA 

A 


Shear area about x-axis 
Shear area about y-axis 


= I. 


PP 


/ y^ dA 


= I„„= J dA 
qq ^ 


Mixed shear area = I 


■rp“ 

= 

pr 

s 

A 

py dA 

— 

= 

/ 

p X dA 

Tq 

qr 

A 

Jr 


= = 

J 

X y dA 

■pq 

qp 

A 



Polar moment of inertia of the cross section about the 
center of shear = I, 




/ (p^+q^) dA 


(13) 


=3 
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Using Ecins . (9) to (13), Eqn. (8) becomes. 


T-U= S 


I i A * Ig 9: 


yB“x«z ]-E [/' 


* 2 Ipq(<-ey){u;,. 9^) 

* 2Irp(4"ex) (<-«*)" 

+ 5J]]ci2 (1^) , 


2,1 .3 Basic Equations : 

Applying Hamilton* s principle to Eqn. (14) and taking 
the variation one gets the seven equations of motion 
associated, with bu^^, 6u^, 6©^, 6Uy, 6©^, 6©^, 60 
and the corresponding boundary conditions, as 

(EA u') - If (pA u^) ■ =0 (15) 

IqqC“x-®yj>^ k<-^ 

h y^k) 

= 0 (16) 
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( 17 ) 


- h^P^ XgO^) 


= 0 


( 18 ) 


fe(^Ixx®x)-OSp‘4**x)-°IpqK- ey)-°Ipr(®z-«'>-|¥<Plxx®x) = ® 


(19) 


= 0 


|^(GI^^(e' -2f))+|^(GI^p(Uy+©^))+|^(GI^q(u^- 9^))+ |7(GI^<Z5) 

- ■k'^P^B^z)"' It^P'^ k^P" ^B^x^ 

.... ( 20 ) 

|^(E 0 )^ GI^^(9;-0)+ GI^p(u^-. ©^)+GI^q(u;- 9y)- 01^(9;- 0) 


k (Pl<. ««) 


= 0 


.... ( 21 ) 

and boundary conditions as 

SA 6ug^ |1 = 0 (22) 

t G Iqq(<-®y)-^ ^ ^ ^ 


El 9,' 69, 

yy y y 


= 0 


.... (23) 
.... (24) 
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n 

KJt 




«y) ^ 




1 


0 


0 


El 6 q 

XX X X 


= 0 


.... (25) 

.... ( 26 ) 


[G + G I^„(u;+0^)4G ©„)+GI^ (25] 6©, 


rp' y X 


rq 


0 


El (25'^6Qi 
ww 


= 0 


.... (27) 
(28) 


In the case of open profiles Saint Venant torsional constant is 
to be used for in Eqns (20), (21) and (28). 

In the above equations' 

EAu^ = Extensional force 

3 . 


[ GIq^(u^-Qy)+GIqp(Uy+9^)+GIq^(©^- 0 ) ]= Ghear force in x-directio 
Elyy 9y = B G ud i n g moment about y-axis 

[GIpp(Uy+©^)+GIp^(u^- 9y)+GIp^(92-(25) ] =Shear force in y-direction 
El ©^ = Bending moment about x-axis 


[GI^^(©' -(Z5)+GI^p(Uy+©^l^^(u'^- 9^)+.G1^«5 ]='..-arping twisting 

moment , 

EI^ 0 ' - Bimoment 

2.1.^ Non -Dimens ionali sat ion of equations and B,C. 

The following parameters (denoted by capital letters) 
are used to non-dimensionalise the equations of motion, 

Eqns. (15) to (21) 





Here the quantities with subscript o represent the section 
properties at 2 = o (one extreme end of the beam) and the 
quantities with subscript n are the non-<iimensionalised 
section properties , 

V/ith these parameters equations of motion become 


.... ( 30 ) 



• • • « 


( 31 ) 



1 


Se?_ Se^_ 


qq 


Se- ^Pn y 

qp 


se^ ^ ^yy- 'y 


♦ * 

©., 


= 0 


qr 


.... (32) 


Rp'^ Rp'^ 

^L. (I (U^+9 ))+ '— (1 (U^ - 0 )) 

pp pq 


Re 


XX 


s© 


F' 

pr 


I- (Ipr„(<-^>))- A^VXgA^^ 


= 0 
.... (33) 


az^^xx^ ®x^~ Se^ Se^ 

PP pq 


Se 


h ^r.-r (%- § )-P, 


2- -pr^-z- T 


pr 


= 0 
(34) 


Re 


■&)w a 


rrr ^ fe ^^rp„^VSx» 


Se 


rr 


rp 


n 


Re^ 


. _:»l_ (I („;.«)). ^|_(i^$ ) 


Se 


2 az rq 
rq 


n 


(£> e 


n 


- \ ®z-*- ^ ^x 


= 0 
(35) 
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^ - 1 - I 

^®rr ^ ^ Sef_ ^ ^ 


ri- W (<- Sy)- A Id (<- ^'I-^J 

i>e ^n ^e n n 

rq 


^here 


_2 . Se^ 

2 » ^®qq 


1‘ 


Co2 

be 

qp 


GI^r^ 1 


P 0^2 

"Tl * ®"qr= 


GI 1' 
qr-Q 


; se: 


GIr^r, 1“ 
PP 

o 

El 

XX _ 




1' 




GI 1“ 




2 * 


GI..^ 1’ 


A T 2 


.... (36) 


GI 1‘ 


GI^^ 1' 


.... (37) 



Similarly using the Eqns . (29), the boundary conditions 
(£qns, (22)to (28)) become 




= 0 (38) 




+ -4— (<- i ’ ) ] ^u„ 


oe 


qr 


qr^' z 


®y 


= 0 ....(39) 


= 0 ....(AO) 


'■ Se2 3^2 


PP 


pq 


^ Ir,,. <-<-P ) ] 


Se 


I 0* 60 

XX X X 

n 


pr 

1 

o 


■pr ^ z 
^ n 


= 0 ....(41) 

= 0 (42) 




rr 


rp 




fi' 4§ 

' n 


Se^„ ■'^n 
rq 


<»)e n 


0 


(43) 

=* 0 ,...(4-4) 



For a prismatic beam (uniform in cross section) the quantities 
v/ith subscript n will all be equal to unity . Equation of 
axial motion (3o) is independent of other equations and its 
solution is well known, hero onwrsrd this equation is not 
considered in the analysis. 


2,1.5 FEM Equations : 


Equations (31) to (36) are the six coupled governing 
differential eauations of motion for thin-walled beams. 
Solution of those coupled equations is obtained by Galerkin 
finite element method in this work. 


Beam is divided into number of elements called finite 
elements, Fig, 3 , The displacements over the typical finite 

element, e , Fig. 3 are approximated as 


U, 


U. 


(e) 

X 

(e) 


©. 


(e) 


X 


© 


(e) 


© 


3> 


(e) 

z 

(e) 


+ a-,Z + 
o 1 

Co + c^Z + 
Cq + c.,Z + 
do + d^Z + 

e^ + e^Z + 


= ILJ 
= lmj 

= L K ] 
= 10J 

,= IPJ 

-IQJ 


Fy] 

l®xl 

ifi 


(ne) 

(ne) 

(ne) 

(ne) 

(ne) 

(ne) 


• • ♦ # 


(45) 







Fig. 3 Finite Element Divison and the Typical Element 
of a Beam. 







Substituting these in Sqns, (31 )ta- (36), six residues 
U © U e 9 $ 

^ . Re y, Re y , Re Re ^ , rI 


are obtained . These 
six residues are minimized by Uin (iali’i-kin method, i.e, 
six residues are multiplied by L^, N^, 0^^', and 

and integrated by parts to reduce the order of derivatives 
wherever possible. For brevity sake, this opor;ition is shown 
here only for Eqn. (31), thus 


Rel. ,Ce) (e) 


(C - e. 


'i .-2 '-qq„'-x -y 


)) 


oe 


qq 


^2 qq ^ X y ^ 

qq 


o o Se^ 


Re 






qp 


qp 




2 

qr 


ih h 


n 


Re, 

Se 


2 

¥ 

qr 


n 


- J h An ’4®’ - Ih '^o An dz = 0 


( 46 ) 

Observing Eqn. (46) and five similar equations, it is 
seen that highest order of derivatives of displacements 

Uy» ^y* $ in the integrals are of first order. 

Thus compatibility of these displacements is needed at the 
nodes of the element. It is also seen that highest order of 
derivatives of these displacements in these six equations' are of 


first order. Thus completeness of these displacements and 
their first derivatives is needed. Hence the compatibility, 
and completeness requirements for all displacements are 
identical , and 

Ll J = [mJ = Ln j = [0 J = [pj = [q J = [nJ (A 7 ) 


The above requirements are satisfied v/ith elements. 

But it has been seen that this element leads to ill- 
conditioning in beams and plates with shear effects, [23,24] 

To avoid this ill -conditioning elements are used in this work. 
For this element, shape functions are 

L N j = [N^ N 2 N 3 J (48) 

where 


N. = 1 - 3 ^ + 2 
' h^ 





N4 = 


T?- 

KT 


(49) 


and the nodal variables are U , U ' and similarly for 
other displacements . 

Equation (46) and five other similar equations are 
put in the matrix form, using Eqn. (45), thus 




V/here 

h h 

[M.)] = J (nUnJ dz ; frig] = J {NRhJ dz ; 


[Mj] = Re^y ; INJLNJ dz 


[M4] = - / A^ X|^ 1 .N 1 LNJ dz,i [M5I . Re^ 5 lN}LNj dz; 


h 


[Mg] = Reg / I3 [N{ LNjdz j [m^] = Reg^ 7 ^NflNj dz 


o n 


2 h 

o n 

2 h 

S 

0 n 


.... ( 51 ) 


and 


Re^ h 


Rei. h 


“^^qq ° ^ 


Se^ o 

qq 


Re„,, h , , 

‘ ^op ^NilNjdz; [K4] 


oe 


Re, 


qp 

2 h 


o ^Pn 


F ^Iqp^NUNjdz; 
0 ^Pn 


Se 


qp 

2 


[%> ^ i VjNllN Jdz ; [Kg] 


-ReL. h 


JU. 


qr 


n 


Se! 


^ ^NILnJ dz ; 


'qr 


^ qr^ 

o ^ n 


Re. 




^h> t¥ ^ ^qq t^lLN J dz -fRe^ / T ^kt}Lk J dz ; 

3^ o o ^ ^ Ti 


’■Re 


qq 

2 


-Re"^ h 


[ K8> i Iqp tN} IN J / W ‘iNl LNJ dz ; 

be o Se„^ o 

qp 


qp 


-Re.l. h 


,^0> i V„ ^NUlijdz ;[Ki^]- ^ 


■3e„„ o n 
qr 


Re^„ h , . 

^ ^ ^nr i^} ^ 

o ^ n 


qr 


31 


[Ki 2 ] = 


^ ? I tN}lN]dz ; [K^j] 
PP 


Re 


XX 


ise 


PP 


/ i„ JLnJ dz ; 


0 PPn 


h / / 

Pl4> > V <M}U»Jdz ; [K15] 

OB o n 
pr 


-Re„„ ’n 


‘"'^nr 


22 s f T J * 

^ -^prn 


Re5„ h o ^ ' r I 

[^163= -r? ' Ipp LN]dz ^ Re^ 1 iHj [tlJdz ,■ 


iie" o 
PP 


h 
r 

0 '"'n 


Re^ h , ^ ► . 

iX > t'Ss] “ Ipr « 

bepr o n oSp^ o n 


R^^„ h 

(^^9> ^ ^rr 

^ Sef„ 0 




(n'} [N Jdz fCKgo] = — s In} lN]dz 


rr 


a>e 0 n 

^®( 0 W ^ 


-irP i 

‘• 5 ev>^ o ^ 
rr 


F2i]= 5 "rrjN} INJ ■i"- ^ \ \ ^*'0 XN] d- 


rr 


" ‘>4 ; (N} LN'J 

o n * 


( 52 ) 



1^1 1 =' 


1 , 


Se^:, 

qq 


be“ 

qp 


Se„^ n 
qr 


®y)* Se2/«’„'''y*®x)-^V„^®x- 

qq qp q^ 


l^'z! = 


‘yy„ "y 1^ 


^yy ®y lij 


{*'3}= 






(e'-$))! 


(f’4l = R®; 


I 0 

■^XX^ X 


I 0 

XXj^ X 


0 
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-V-^rD ^^'+0 )+ -^1^.. (U' ~ %) 

^®rr 2®rp ^ ^®rq ^ 

fi>e n 

0 

^c:«2 ^rr^^®z"^ .^-^-^..,i^p^(Uy^0^)+ ^._2 ®y^ 


Se“ ""n 
rr 


Sef. 


rp “■ ^’^rq. 


+ f ) 

we n 


= Re. 


^'000 9 


toor \ 


n 

0 

'■CiMjJ- 


$ 


n 


0 


'.... (53) 


[uation (50) in a concise form is 


r T (®) r"i r r ^ (ne) .(ne) 

[m] {u} +[k] {u} = {f} ....(54) 


lere 


(e) 

[mG = 

^ ^e) 

uf ^ 





mass mtrix of the element 
stiffness matrix of the element 

displacen^nt vector of the element 
boundary force vector of the element 



Assembling for all elements and applying the boundary conditions, 
equation of motion become 


.. (n) 

[M ] (U} . 


(n) 

[K] iu} =[0j 


.... (55) 


2.2 STATIC rdJALYSIS ; 

2.2.1 Basic Equations : 

The governing equations for the static analysis 
are derived (including shear deformation effects). 

For this case potential to be minimised is 

V = U - Ug .... ( 56 ) 

Where U is the strain energy given by Eqn. (7), u is the 

@ 

work done by external forces and is 

1 ^1 


"*y“y® Vx‘^ (.z-I^J*l-^. S^l>iz {z-z ) 

+ )- 1^,0 (^-^0)) dz. 

v/here 

Qvj’Itt loads per unit length of the beam in 

the x,y and z directions respfectively • 

®t’®x’™y ~ applied torsional moment bending moment 

(about x), bending moment (about y) and . 
bimoments per unit length of the beam. 


35 


Q^,Qy and = concentrated load in the x,y and z directions 

respectively 

Mt,M^,I'Iyand ]\!^= concentrated torsional 

"bending moment (about x), bending moment 
(about y), and bimoments, 

6 (z - z ) = Dirac delta function for concentrated loads 

and moments 

Minimising the potential, Eqn, (56), one gets the 
governing equilibri\im equations as 


l^(EAu^) + ^“^u V ....(58) 

a 

fe(GI^,K-9y))- fe(GIqp(u;* e^)) 

0 ....(59) 

fe(EIyye;).GIgq(4-6y)^GIqp(u;+e^)^GI^^(9;-5() 

+ niy+My. ^(z-Zq ) = 0 ....(60) 

|5-(GIpp(u;+e^))^ |j{GIp^(u;-ey)) +|^(GIpr(9;-8!)) 

+ qy+Qy ^(z-z^ ) =0 ....(61) 




|--(si o')-GI (u^ + 9 )-GI (u" - 0 )-GI (g'-CZI) 
3z'‘ XX x-' PP^ y pq^ X pr^ z 


+ m^+M^6 (z-Zq ) 
IK 


..( 62 ) 


fe(GI^(e; -«!))- fef«rqK-®y^) 




..(63) 


|-(ET 0)+ GI (o'"-(^)+GI (u*+0 )+GI (u'" -9 ) 

3z^ rr'- z rp^ y x^ rq^ x y' 


rr' z 


- 91^(9^ -5^)- m^-f^6(2-z^) 


= 0 ...,(64) 


2.2.2 Hondimensionalisation : 

The governing equations (5s)fco(64) are nondimension- 
alised, using the parameters of Eqn. (29), and are 


Tf (^v>U„)+ q^ + Q., 
oz ' n a'^ ^z z 


= 0 ....(65) 


1 3 


1 3 




* + Q; 


0 ....(66) 


Irdvv -V- ^nti 

^ YYn y Se^ ^ ^ Se*^ ^ ^ 

99 9P 


+ — ^ Iqr * V^V 


= 0 ....(67) 
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^ Irdr.. 9y)) 


Sg2 32' PPn y 

PP 


oe 


pq 


+ — i- It (©J- 

S 2 3z P^J^ Z 
®pr 


))+ Qy + Qy = 0 


® <)- -V^nn (“v-"*,.)- (U^- e„) 


Jz'-xx^-X' se^„'PPn''y ■='' Se^ >■ 

^ PP pq 


-^pr 
'pr 


V ^®2' f 5 “X * % 


= 0 ... 


Se 


2 az ^ Se^ ^ 


rr 


rp 


n 


1 a 


fe(Irq„<«X-9y))- fe (\ 


= 0 


rq 


(9;- §)- ^Irp (V®x) 


'n "n 


Se" '^n 
rp 


+ Sy)- ^^2 ^<9^- J)-ma - H 


Se"^ 

rq 


0 ... 


Whf^rp Se^ . Se^ , Se^ , Se^ « S^ , Se^ . Se^ Sp2 
wnere ^p, q^,, «pp* , oe^^ , oe^ 


and cs'e^ are defined by -Eqn, (37) and 


.(68) 


.(69) 


(70) 


.( 71 ) 
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'V 






'^x “ 


0x1 ) 


El, 


yy. 


m,, = 


V' 

El » 
yy^ 


My 


Myl‘^6(z-ZQ ) 

y~ 


yy. 


% 


ET~ » 

xx-o 


Q. 


Oyl 


El 


XX, 


mx 


m~ 


XX, 


% 


6(z-Zq^) 


El 


XX. 


m,|,l 


“I " Ei;^ ' 


M^l '^(z-Zq ) 


o 


«1q = 


n^l 


I M ~ 


-SCz-z^) 


(72) 


2 .2 . 3 FEM Equations ; 

Equations (66) to (71) are the six coupled governing 
differential equations for the static case. Solution of 
these coupled equations is obtained by Galerkin finite 
element method in this work. 
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The displacements over the typical element, B^ig.? , 


are aonroximated as in (^5). Substituting Eqns.(45) 

' U Q U 

in Eqns, (66) to (71), six residues Re , Re ^ , Re 

6„ © d 

Re , Re and Re'^ are obtained. These six residues are 
minimized by the Galerkin method. For brevity sake 
this is shov/n here only for Eqn, (66), thus 

/(e) (e) h h 

))1 

o 


N -y- (I (U -© )) 

qq. 


1 / ,(e) (e) 

- ^ ^^x - ®v 

o se'^. ^ q^n ^ y 


qq 


1 , , 

N. — 4-(I (U + 

1 qPn^ y 

qp 


9. 


x 


)) 


h 


.(e) _ (e),, 1 “ . y ^ - 


o oe 


qp 


,(e) (e) 

\ 

X 


, (e) (e) 


Se ' ^ n 

qr 


- I 


i q^n ^ 


0 0 oe 


qr 


h 

/ q^ 


h 

dz + Qx 

o 


0 


. (73) 


By observing Eqn. (73) and five similar equation for 
compatibility and completeness, one can see that elements 
can satisfy the requirements. But it has been seen that this 
element leads to ill -conditioning in beams and plates with 
shear effects, [23,, 24 ] . 'I'o avoid this ill conditioning 

element are used in this work. 6o the interpolation 
functions of. Eqn, (48) are used. The modal variables are 
and similarly for other displacements , 



Equation (7j5) and live othor equations 


are put in the matrix form usin^ Eqn. (45), 


T (e) .(ne) . (ne) Ane) 

[k] [uj = {q\ + [f\ 


.... (74) 


r 

Where [ K ] = stiffness matrix of the element 

. .(ne) 

\_U ^ = displacement vector of the element 

f v(ne) 

“iQj = Load vector of the element 
f Une) 

\FS = Boundary force matrix of the element 

r ie) . (ne) 

Matrices LK J and are obtained from Eqn, (50), 

by substituting Re^ = Re^y = Re^^^ = 1 . 

(ne) 

And load matrix {q| is 


0 

7 iOy dz + / My dz 

7 {n| qy dz + ^ {N^Qy dz 

o 

h ^ r I 

; {N} ray dz + / (WfMy dz 
o o 

/ {n} mm dz + / (n Mm dz 

o o 

h . . 


(ne) 



1 

( 


1 


] 

rt 

iQijl 


^5! 


(“6} 

- 


(75) 
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Assembling for all elements and applying the boundary 
conditions, the governing equations for static case become, 

[k] {u] = .... (76) 

2.3 BUCKI^inG ANALYSIS : 

2.3.1 Basic Equations ; 

In this section coupled flexural -torsional buckling 
equations of a thin walled column are derived from variational 
principles, These equations are classical equations and have 
been derived using equilibrium method by many researches 
See [17,2 ] . 

The inplane displacement, v and w of a point A on the 
cross-section of a thin-walled column. Fig. A-2, are given 
by Eqn, (A-9) . 

Using Eqns . (A-3) to (A-6), iiqns.(A-9) become 


V 


w 




(77) 


where 

u , u = the displacements of pole B in the x and y-directions 

^ ”b 

respectively 

Xg,yg = coordinates of the pole B. 
a = angle of tangent at A with the x-axiS' 
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= the perpendicular distance from the point 0 to 
the targent at point A 
*= (x^,^ sin a-y^ cos a) 

= The perpendicular distance from the point 0 to 
the normal at point A . 

= cos a+ sin a) 

“ Coordinates of the point A 

The out of plane displacement of point A is given by 

Sqn, (A-22) . Neglecting the effect of shear strain, it becomes. 


u 


^a- 


/ 



/ 





(78) 


Here the subscript b has been dropped. 

The point 0 is taken at the centroid of the cross section 
and the point B as the shear center in this v;ork, Fig, (A-2) . 
Here onwards subscripts B of the displacements and their 
derivatives of the shear centre B are dropped. 


Normal strain is 


e 

zz 




The contribution of 


1 

1 



2 

) to the e _ 
' zz 


is negligible 


compared to others, thus 


e = 
zz 




( 79 ) 
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VJith this strain energy becomes 
U = ^ / (E dv 


V 


ZS ' 


. ♦ / TT r ^ / / 3v\ 2 / BViT \ 2 \ / cju ov \ 2l j ■ j 

; ; |E [^ + ^((^) +(^) ) ] jdAdz 


o A 


•*.. (80) 


Substituting the Eqns . (77) and (78) in Eqn. (bO) and 
neglecting the terms having highor orcior go'-'ors of the 
derivatives one gets , 


U = 


1 


1 

/[ 

0 






“4 




+ Au; 2yg AU^ Auj, S^u^) 

+ GI^ ] dz .... (81) 


Terms Ig and are defined in Eqn. (13). 
Now treating SAu^ as a constant, axial compressive force, 
-p , Eqn, (81) becomes 


U = 


1 

2 


1 


/ fe (Au ^ ^+I ^+I u" ^+I 

'■ a yy x xx y 


0)W 



- zpyb u;^e'^+ apx^ e- + gi^ ] dz 



li # * • 


( 82 ) 
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Minimizing this potential U, one gets the four buckling 
equations associated with du , du , du , dQ as 


up 




2 
3z 

3^ 


32 


2 3z 32 ®2^ 


h ^"TT ®z) 

3z 


and the boundary conditions 

il 

= 0 ; 


EAu„ 

6u' 

a 

a 


u" 

yy 

X 


= 0 «... ( 83 ) 

= 0 .... ( 84 ) 

= 0 .... ( 85 ) 


*" Tz ■*' 3z ^x^ ” ^ •••* 


= 0 * 


... ( 87 ) 


= 0 


.... ( 88 ) 


SI u''du 
XX y y 




= 0 


* 4 » # « 


( 89 ) 





tf i 


El 9 60 

2 Z 


Pli 




-oxg u^- Gi^eJ) 6e^ 


= 0 


(90) 


2,3,2 Noii-dimensionallsation : 


The buckling equations --‘jns (83)tD(&5) are non- 
dimensionalised, using the parameters o.f ann. (29) 


These become 



(91) 


= 0 .,,.( 92 ) 


Hy fia ®z) 

^ ■ §2 


0 • • . • (93 ) 
0 (94) 


p R&^ I 

( I ©”)- — q')+ Re^ ^ (P ©*) 
c)Z n we oZ n n 


- 52 up ' ° 


V/here ReJ^, Re,^„, Re,?,.' “e^ and Re| are defined in 

XX' yy wo) 45 

Eqns . (37). 
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The boundary conditions are 

= 0 ,• 


“a 


. ; . .(96) 


I U 6U 

YYn ^ 






= 0 


o 


tt „ / 

I U 6U 

y y 


I Q <5© 

(WQ^ Z 2 


(97) 


0 ; (Re2^'-yi^^ ^Vfu'-PXj, 


= 0 


n 


PIb 

0 i (P®(^co az^^COO) "T 

n ^n 


0 

(98) 




PYr,U. -PXwU.- 




<«)e n 


1 

0 


(99) 


2,3.3 FHM Equations : 

Equations (93) to (95) are the three coupled governing 
differential equation of buckling of a thin v/alled column. 
Solution of these coupled equations is obtained by the 
Galerkin finite element method . 

The displacements U , U and 0_ over the typical 

X y z 

element, Fig. 3, are approximated as in £qn. (45), 

Substituting Eqns. (45) in Sqns . (93) to (95), three 
u u e 

residues Re , Re and Re are obtained. These three 
residues are minimized by the Galerkin method. For brevity 
sake this is shown here only for Eqn. (93) thus. 
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'»(®) 1^ p t 

-Re"^ N, I U 

o yy i yyn 


^ h ,, (e) ,(e) ,h h 

+ Re£. / N, I.... U„ dz + N, PU 

J- 


yy o ^ yyn 


,(e) 


- J N, PU dz 


, / \ .h h - , 

+ N. PY^ - ; N. PYt 5 dz = 0 

1 B z 1 ^ ^ 1 B z 


(e) 


... (100) 


By observing Eqn. (100) and two similnr cuu.-.tions for 

compatibility and completeness, one can see that elements 

for U„, U . and 0 satisfy the requirements. The interpolation 
X y z 

functions of Eqn. (4s) are used. The nodal variables are 
^x’ ^y» ^y » ®z node. 


Equations (100) and tv/o other similar equation 
are put in the matrix form using Eqn. (45), 


[K^] [O] [O] 


[0] [K,] [0] 


[0] [0] [Kj] 


{Ux\ 


(ne) 


(.'^yl 


(ne) 


l®z] 


(ne) 


Fl^;] [0] [Lj]' 


}-F [0] [q ] [ Lj] 






(ne) 


A^yl ’ / 




(ne') 






{^31 


( 101 ) 
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V/here 


.2 H T im'i dz . [k^]= Re^ / 


'yyn 


[Kj] - ^ LNJdz + {n} Ln Jdz 


[m = ^ l^nU'^J dz ; [L2> ^ 

O ^ O 


....(102) 


i 

[L3]= - / Xg l_Njdz i [L^] = Re| J {.Wj UJ dz 


and 


(^i} = { 


(ReB.. I7 (I,„, UJ + PU^P Y„s') 


■yy ^ ' yyn ='" 


Re^ I 

\yn “x 




-CP^yy h ^hvS^* ^ 


- ^y V, < 


(p?V 


(Re^ Iit(I U ) + P U'- PX^e^) 


1 




-Re^ I 


n 

// 


U 


XX xXj^ y 13 

■^^®xx I7 ^^xx^ Uy)+ P Uy-P Xg 0^') 

- Re^ I u' I 
XX xx^ y ly. 


I (103) 
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^''3) = 




Re,^.. o IB 

T Q* -f. Re^ -r— <2 -n q' 

COG n n 


p Xg U- * P Yg u; ) 


- 


? // , 

Re I Q 
coco z 


( Iz'Ioc ®z)* % ®z- P ®z 

n (oe n n 




. o Bv, 

ww T r,' t>_ 2 n. 


+ P Xv. U, 

3 y 


- P > 


Re^ I 9" 

ww z 


Assembling for all elements and applying the boundary 
conditions, Eqn. (101) is solved for the critical buckling 
loads . 



2.4 PLANE FRAME : 


In this work plane frames (lying in zx plane) 
subjected to transverse loads (y -directions) are considered, 
plane frames are formed by joining thin walled beams . FEM 
equations for a typical element are first obtained in the 
local coordinates (z x) and then in the global coordinates 
(zx) using a transformation matrix.Shear effects are 
neglected in the frame analysis . Frames are studied for static 
and free vibration cases. 


2.4.1 Static analysis : 

2.4.1 .1 Basic Equations : 


As the plane frames subjected to transverse loading 
are under study, the governing equations for the beams of 
these frames are Eqns . (61) to (64). Neglecting shear, these 
four second order differential equations are combined to two 
fourth order differential equations, 

3 m, 




El — - 7 ' ^ = q +5 6(z-zu ) + 

xx|pr S y ^ 3 - 

4— 2 — 

3 3‘^», - 3m,. 

r:3r - 7::r = •■■■ Co5) 

3Z 32 Z 32 


X 


.... (104) 


The over dash refers to the displacements and 
rotations in the local axes. 

The above equations are valid when the loading is at 
the shear center of the beams. When the loading is at the 
centroid of the cross section of the beam, the equations (104) 
and (ip5) become [2 ] . 



X. No. 


A: 


TTPTZZJ 
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0 *0 


4— — 

3 u 3m 

El = q + S 6(z- 2 )+ — i 

XX ^ y ^ 3£ 


(106) 




3^© 


3m. 


EIou -Old ^ — -5x(yq-yB> 


5y(Xq-Xg) 


.... (107) 


where (x ,y ) the coordinates of the point of loading and 

M. Vi 

(Xg, Yg) are the coordinates of the shear center. 


2,4,1 .2 FEM Equations : 

Equations (106) and (107) are the two differential 
equations for the static analysis of a beam of a plane frame. 


The displacements over the typical finite element, e , 
Fig , 3 are approximated as 


(108) 


©, 




Substituting these in Eqns. (106) and (107), two 
residues Re , Re are obtained . These two residues are 
minimized by the Galerkin method. 

Study of compatibility and completeness requirements 

- (e) 

shows that elements are needed for displacements Uy^ ' and 
thus 


LN J = LM J= IN J 


• # • « 


(109) 
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( 113 ) 





-El 5 
(i)W z 


(-EIoh^Sz Id s;) 


EI,„, 5, 

ww z 



where 

mm ^ 

^y " Bending moment ; - EI^ Uy = Shear force. 


— << ftf / 

"^^oxo^z Bimoment ; ^d^ Warping Twisting 

moment +St Venant 
tisting moment 


= Total twisting momeni 



Equation (110) in a more concise form is 

= 1S( [P] .... (114) 

where 

I" K = Stiffness matrix of the element in the 

local axes . 

[uj Displacement vector of the element in the 

local axes. 

[ <3j Load vector of the element in the local axes. 

1 ?'} Boundary force vector of the element in the 

local axes. 

2,4.1 ,3 Transformation to the Global axes . 

The displacements and their daxivatives of a beam 
element of the plane frame in the local axes (x,z). Fig, 4, 
are transformed to the global zes (xz) as [18,6 ] 
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Fig. 4 Local and Global Axes of a Beam Element. 
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or {ul = [t] .... (115) 

where a is the angle between the element jk and global axis z . 

Matrix displacement matrix in local axes, 

[T] = transformation matrix, 

displacement matrix in global axes. 

Similarly forces in local axes are related to 
forces in global axes, 

=[T] ....(116) 

Using the Eqns . (115) and (116) , Equation (114) becomes, 

.... (117) 

Where 

[ K]^®^ = [ K ^®^ [^1= Stiffness matrix of the element 

in the global axes, 

£Qj load matrix of the element in global axes, 

/ \ 

[ F[ '^= boundary force matrix of the element in global ax< 

After assembling for all the elements, using the Eqn. (117) 
and applying the following boundary conditions. 

Pined end Uy= 0 ; EI^ u"= 0 ; ° ^ ° ...*(1^8) 

Fixed end u.y,= 0;Uy = 0 ;©2 = 0;©2 = 0 ...,(119) 

Free end El u''= 0 ; El u"= 0 : El ©"= 0 ; 

(ET ©''-GIj ©^ )= 0 
(*) 0 ) z d z' 


( 120 ) 



one gets. 


[ r iu(" = wi" 


( 121 ) 


Equation (121) is solved for the displacements in the 
global axes . 

2, 4. 1.4 Stresses in the plane frame . 

After solving the Eqn. (121) for displacements, 
one can find out the nodal forces and moments 
from the Eqn. (117) as 

.... ( 122 ) 


Using these moments and forces, the total stresses 
(Normal and shear) at the nodes are calculated as [2,5] 


M. 

Total normal stress = a = >»■ 

z I 


XX 


y o(sg) 


(123) 


XX 


( 0(0 


Total shear stress 


Q I 

- - V X 

^zs"* t X 

XX 




0X0 


.... (124) 


where 

^xx “ Bending moment = E 

M,. * Bimoment = -E I 6* 

(0 (0(0 z 

(o(Sg)= Warping function, see Appendix III and c 5 ] 

I 

Q„ = Shear force = - E I u*! 
y y 
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I = Statical moment about x-axis = <1 y 

^ A 

t = Thiclaiess of the profile 

Mj. = Warping twisting moment or flexural twist^-EI^e^ 

I (s^)= Statical sectorial moment = I w(s„)cLA, see 

6) £1 Jh 

Appendix III and _ 5 ] 
= Saint Venant twisting moment = GI^ 

= Saint venant torsional constant of the section 


2,4,2 Free Vibrations : 


2. 4, 2.1 Basic Equations : 


Out-of -plane vibrations of the plane frame are 
being studied . The governing equations of motion for this 
case are Eqns . (18) to (21). These fovir second order 
differential equations can be combined to two fourth order 
differential equations , Neglecting the shear one gets these 
equations as, 


' aV 

El — yX _px -L_ 

^■^xx ^-4 ^-^xx 

9z 3z at 


♦ PA -pAXg 


a ^5 

z 


= 0 


« ♦ • • 


(125) 
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aV 


a^5„ 




-Id -PI.. ^pIb ^ - p^^ ^ 


( 126 ) 


2. 4 , 2. 2 FEI-1 Equations : 

Equations (125) and (126) are the two differential 
equations for the free vibration analysis of a plane frame. 
These two differential equations are solved by Galerkin finite 
element method . 

The displacements over the typical finite element, 

e , Fig. 3, are approximated as in Eqn, (108). 

Substituting these in Eqns. (125) and (126), two 
u 5 

residues Re ^ , Re ^ are obtained. These two residues are 
minimized by the Galerkin method. 

Here also C.j elements are needed to satisfy the 
compatibility and completeness requirements. The shape function 
for the elements are as in Eqn. (49). The nodal variables 
are Uy , u^ and 

The minimized residual equations are put in a matrix form, 
using Eqn, (108) as, 

; M](®> tS! - iPi 


• t • # 


(127) 
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Matrices [K [uj {Pi are same as in 

Eqns , (111) to (114), and matrix [M is 



Mass matrix of the element 
in the local axes 




r Mz 


(e) 


J 


(128) 


where 


](®)= pA ; {N? 1 NJdz +pl^ f In] [N'jdz ■ 


(e) h 

- M^l = -pAXg / {n1 _ n: d£ 

o 


(129) 


li In 

r M^-j (®)= pig / {Nj i N! d£+p / iN'n N'jdf 


2. 4. 2, 3 Transformation to Global axes : 


Equation (128) is applicable to the elements in 
the local axes (x,z),Blg.4* This equation can be transformed 
to the global axes (x,z), following the procedure of 
section 2.4.1 ,3, as 

(e) (ne) 


[M] {u] +[K1 


(e) (ne) (ne) 

lu] = [Fj .... (130) 
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(e) (ne) ,(ne) 

Here [K] , ^uj and iFi are same as in Eqn. (117) and 

(e) J - (e) 

[Mj “[T] I’M] [T] = Mass matrix of the element in 

the global axes, and 

[ T ] is the transformation matrix, and is same as in 
Eqn. (115) 

After assembling for all the elements and applying 
the boundary conditions, Eqns, (118) to (120), one gets, 

[m] [u? ^ + [K] lu]” = [O] .... (131) 

which are solved for the natural frequencies of vibration. 
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CilMpT^Lit- III 

RSSULTS AND DISCUJCIO:J 

This chapter deals with the niimerical results, 
obtained by Galerkin Tinice elemtsnt mrthod, of the coupled 
equations of thin walled beams for free vibr; tion, static 
and buckling cases, for various boundary conditions. The 
results foi' free vibration and static analysis of a plane 
frame are also presented. 

3.1.- OF FIl^gTE ELEMBNTS : 

First of all one should decide the number of finite 
elements to be used in the analysis . For various boundary 
conditions, the problems v/ere solved for two, four and six 

finite elements . It was observed that v;ith six finite 

elements, first five frequencies v/erc feirly accurate in all 

the cases. Sample results foi’ pinned -pinned thin-walled beam 

of unequal channel cross section, Fig, 7 , are shown in 

Table 1 . 

3.2 FREE VIBRATIOMS : 

For a beam of* doubly -symmetric cross-section, shear 
center B and the centroid 0 coincide. 

Thus 

X,. = Yt, = I'=I sI =0 ....(132) 

^ B pq qr “ pr , " 
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Table 1 : Comp 'risioii of >''Xact and ‘fl solutions for 
a pinned -pinned boom of une''ual channel 
Section. 


Mode No . 

n 


Natural frequencies 


Finite element method 

sol u Lion 

Nxact 

• solution 

[ 9 ] 

2 elements 

4 elements 

u elements 

1 

400.10 

398.93 

398. r.S 

398.84 

2 

592.73 

590.60 

590.47 

590.44 

5 

161 5.44 

1469.26 

1464.99 

1463.84 

4 

1805. 11 

1798.51 

1798.12 

1796.06 

5 

2603.11 

2355.36 

2348.03 

2346.14 
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With thio governing equntioni^ , (31) to (36,-, become 
uncoupled , 

A beam having one axis of symmetry, rhe shear centre lie 
on the axis of symmetry. In this case the bonding vibration 
equations about the axis of symmetry become uncoupled. If the 
x-axis is the axis of symmetry. Fig. 3, then 



.... (^ 33 ) 


Bending vibration equations about y-axis (31) ^nd (32), become 
uncoupled ; and Sqns . (33) to (36) remain coupled. This is 
ioiovm as double coupling of bending vibration about x-axis and 
the toi'sional vibration. 

In the case of beam v/ith no axis of symmetry, shear 
centre does not lie on any axis . Governing equations remain 
coupled and this is knovm as triple coupling of two bending 
vibration about the x- and y-axes and torsional vibration 
about z -axis . 

3.2.1 Mono-axis symmetric beam : 

For this, thin walled be-'’.ra of seiai -circular cross- 
sections having x-axis symmetry Fig, 5, case 1 of [14] 
and case 2 of [ 16 ] are chosen in this work. The dimensions, 
section properties and other parameters are given in Fig. 5, 



y 



Case. 1 

a = 24.5xl6m , t= 4.0x10 m 

E = 68.9Gpa. G=26.5Gpa. f= 2711.039 Kg/ 

Ixx = 92.6x16®m^ Ljtl.52 xl5’*m * 

lB=0.1848xlC)®m^, A =308x10 I = 0.820 m’ 

X =15.5 10 m . Yb= 0.0 m , I pp = 153.938 x 10 

Ipr= LOxlO^^m^, Irr =0.035 xl6®m* Id = 0.8210 xfo^m^ 

Case . 2 .3 

a = 50,8x10 m. t =15875x10 m 

Ixx=0.3269 X 16 m*. = 20.0 x 10 m ® 

Ib= 0.6538 xl6®m\ A = 253.35 x 10 V, I = 0.90 m 
Xb = 0.03 234 m , Yb = 0.0 m , Ipp =126.677x10^2 
Ipr= l,0xl0'm^ Irr = 0.1238x10 Id = 0.2128 x» m* 


! I 



Fig, 5 Thin-Walled Semicircular Section. 
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The governing equations ore (33) to (36) and (133) . 
These equations are solved for pinned -pinned, fixed -fixed 
and contilever cases and the results are given in Tables 2,3, 
and 4 respectively. 

First, case 1 is studied and results are obtained 
neglecting rotary iiiertia and shear. In this case Eqns. (33) 
and (34) become a^. fourth order differential equation, and 
Eqns. (35) and (36) become another fourth order differential 
equation. These tv/o fourth order coupled equations are solved 
using elements with boundary conditions as, 

Pined end = 0 ^ ^ ° 

Fixed end U,, == 0 ; u' = 0 ; =•• U ; 9' = 0 (13^) 

y y z; 

Free end u"' = 0 ; u7= 0 ; g'J = o ; ©1"= 0 

y y c 

Next, results are calculated including rotary inertia 
terras (but neglecting shear). In this case also one gets two 
coupled fourth order equations and are solved using same 
elements and boundary conditions, Fqns , (''34). 

Lastly , results are calculated including both the 
rotary inertia and shear terms. In this case there are 'four 
coupled second order equations. These are solved using 
elements with boundary conditions, 
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Table 3 : Coupled frequencies of a thin walled fixed -fixed 
beam of semicircular cross section. 
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Table 4 : Coupled frequencies of a Thin-walled ccintilever 
beam of semicircular cross section. 
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Pinned end 

c 

II 

o 

; 9^ = 0 

, J ©2 = 

0 f 

o 

11 

Fixed end 

U = 0 

y 

, = 0 


0 , 

1 = 0 

Free end 

^x “ 0 

uVq 

Y X 

-1 * 

n 
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pp 

Se2 

pr 

— ■ U , 



1' = 0 
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U'-i-9 
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(135) 


i^tudy of the results in •'•Voles 2 to A shov/s that 
frequencies decrease due to rotary inertia and further 
decrease on inclusion of shear, as expected. For pinned-pinne( 
beams without rotary inertia and shear, exact results of [9, 1C 
are also given in Table 2, Present FSiM results compare v;ell. 
Error in FEM results increases v/ith increase in mode number, 
as expected. For cantilevorbeam, Friberg [1A ], obtained the 
frequencies including rotary inertia (but neglecting shear) 
by dynamic stiffness matrix method. The results are shown 
in Table A and present FEM results match well with them. 

Me xt^ semi-circular cross-sectional beam, Fig. 5 , 

, case 2j is studied. The ratio of frequencies (^+s/^j{) Tor 
this case is shown in Fig, 6 as a function of 
the frequency of vibration including the effect of rotary 
inertia and shear, is the frequency oi vibration excluding 
the effect of rotary inertia and shear. 




(b) Fixed -Fixed Beam 
Frequencies of Semi-drcular Cross-section 
Beam as Functions of Rexx* 
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Tso [16 ] obtained the results analytically.- Present 
results match v/ell . It is not possible to see any difference 
in the results in this figure, 

Por case 1 , it is seen that decrease in frequency is 
not alv/ays larger for the higher modes than for the lov/er 
modes. This v/as obser-ved by [ 16] also. Ptudy of mode shapes 
for pinned-pinned bean shov/s that first, third, fourth and 
fifth are torsion dominated modes ; and second and sixth are 
bending dominant. Parameters for this section are 


= .000447 , Se^^ = .002326 , Se^^ = 2.9361 
= .00001223, Se^^ = .00016789, Se^^ = .0000716762 


.... ( 136 ) 

fferameters Re„„, Se_, Se_^ of bending are larger than the 
parameters Se^^,, Se^^ of torsion, hence the decrease in 

frequency is larger in bending dominant mode than in torsion 
dominant modes , 

For fixed-fixed beam first, second, fourth and fifth 
are torsion deminant modes and third is the bending -dominant 
mode, thus more decrease in frequency in third mode. For , 
cantilever beam, first, third and fourth are torsion dominant 
inodes, second and fifth modes are bending dominant modes. In 
both these beams, decrease in frequencies is large in bending 
dominant modes than in torsion dominant modes. 




Es68.9Gpa, Gs26.5Gpa. f = 2711.039 Kg/ m’ 

A= 0.01483&62m^ Iwu>= 1.8403 X 10* m®. Ixx =545.263 xlOm* 
Iyy=73-673 X Ip = 618.936x10*^ Ib = 886.91 x 10® m* 

Xb = -0.121793 m . Yr = 0.0567944 m , Iqq = 9X)3224xl6® 
Iqp=l6^*m^ Iqr = 0*283004xl6Vjpp = 5B0644 x lo’ 

Ipr =0.36768 xlO^m! Irr = 0.465625 x l6Vld=225459xl6®m^ 


Fig.7 Thin-Walled Beam of Ur^quKil Channel SecttCMO. 
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3,2.2 Unsymmetric section beam : 

For this, a thin walled beam of unequal channel 
cross section, Fig. 7, of [5 ] is chosen. The dimensions, 
section properties and other parameters are also given in 
the figure. 

The governing equations are (31) to (36). Results for 
pinned -piiined, fixed -fixed and cantilever beams are obtained 
and are given in Tables 5,6 and 7 respectively. 

Firstly^ results are obtained neglecting the effects 
of rotary inertia and shear. In this case the six coupled 
second order differential equations (31) to (36) reduce to 
three coupled fourth order differential equations . These three 

t 

coupled differential equations are solved by Galerkin finite 
element method using elements with the boundary conditions 
as, 

Pinned end = 0 ; 0 Uy= 0 ; Uy= 0 ; = 0 ; 9 ^- 0 

Fixed end = 0 ; ° ^ 0 ; Uy= 0 ; ©^ = 0 ; ©J= 0 

Free end ® ® ^ ^ ^ ^y~ ^ ^ » ^z" ^ 

.... ( 137 ) 

Next, results are obtained by including the effect of rotary 
inertia (neglecting shear) . In this case also there are 
three coupled fourth order differential equations . These 

elements and boundary conditions. 


are solved using the same 



Table 5 : Coupled frequencies of a Thin-walled pinned -pinned 
beam of unequal channel section. 



Table 6 : Coupled frequencies of a thin-v/alled fixed -fixed 
beans of unequal ' channel section 
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Table 7 ; Coupled frequencies of a Thin-walled cantilever 
beam of unequal channel section 
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i’inullyj, ;iru obt-iinud inclufJin,-_ tlie efTecto o£ rot-iry 

inertia, and she‘'.r , In this the governing equations are six 
second order coupled differ jntiul equ;d:ioii 3 . These six 
equations, riqns . (30/to(34), are solved by FEM using elements 
with the boundary conditions as 

Pinned end 0 ; 9' = 0 ; U = 0 ; 6' = 0 ; 0 = 0 ; J' = 0 


Fixed end 


Free end 


^ » ®y =*0 ; Uy = d ; W,,-. 0 ; 02= 0 ; f 


0 


0=0 

y 


’•^x"®y , ^ 


* ? 

qq 


= 0 


o6 


■qp 


■ jQ 


qr 


U/+0. 
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o' = 0 • .f 

X ^ ^ .,,2 .2 

pp pq 


<- i) 


= 0 and .... 0 58 ) 


uo6 


pr 
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0 ; ^ — + 


» r. 2 


rr 


U -f- 0 U - 9 

. , .2£ I 

o o 

rp rq 
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It is clear from the results in Tables 5 to 7 that 
frequencies decrease due to the rotary inertia and decrease 
further with the inclusion of shear as expected. Excluding 
shear ard. rotary inertia exact results can be obtained 
using, Gere [ 9 ] . Fxact results for pinned-pinned case are 

also given in Table 5 and present FEM results match well . 

In the three beams studied effect of shear and rotary inertia 
is least in cantilever and highest in fixed-fixed beam. 
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In n.. - il, in ohcerve^l t: h d.^cr^nse in 

freuin-ncy is not. nl ..nys 1 r£.or for bicVv-T nodes than the 
lower moces . otudy of mode snaues for pinned-pinned beam 
shov/s that second and fifth modes -sre bending; dominant about 
the y-axis direction, fourth mode is bending dominant about 
the x-axis, first and third modes ''re torsion dominant. 
Parameters for this section are 


Re^ = 0.0007943 

; Se^q= 0.003393 

} Se^p= C.3064E+17,'Se^^=0.2707 

Re^ = 0.005879 

; Se2p=0.059 

; 3e^^=0.2268E+l8 ;Se^^=1 .5423 

“ 0*0^^05319 

; Se^ *0.001644 
' rr 

i 3e^p=0.001082 ; Se^q=0.30S33 



(139) 


Derivation of these parameters is given Appendix III. 
Parameters ^®pp» ^®pq ^®pr sending about x-axis 

are larger than the parameters and of 

bending about y-exis. These p-.rameters , Se , Se and 

j y Hni Hir 

^®qr larger than the parameters Se^^ and 

^®rq torsion. Hence the decrease in frequency is larger in 
bending about x-axis dominant modes than in bending about 
y-axis dominant modes. The decrease in frequency of bending 
y-axis dominant modes, inturn, is more than that in torsion 
dominant modes* 
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r 01" fixod -fixed bu?.rii, second, fifth modes are bending 
y-axis dominant modes, fourth mode is bending about x-axis 
dominant mode and firsb, third modes are torsion dominant 
modes, ^'or cantilever beams, second, fifth modes are bending 
about y-axis dominant, third mode is bending about x-axis 
dominant, and first, fourth modes are torsion dominant. For 
both "these beams, decrease in frequencies is more in bending 
about x-axis dominant modes then in bending about y-axis 
dominant maies , Inturn the decrease in frequencies is more 
in bending about y-axis dominant rnocc:-' th‘'n in torsion 
dominent modes. For a bettor understanding of this, the */. 
decrease in freciuencies as a function oC is shown in 

JFig, 8 for pinned -pinned and fixed -fixed beams. Symbol T-j 
stands for torsion dominant j mode, r..-x-j for bending about 
x-axis dominant j-th mode and h-y-j lor bending about y-axis 
dominant o-th mode. 

3.3 STATIC AI-IALYSIS : 

For beams of doubly -symmetric cross-section, Eqn.Cl32) 
holds and the governing Sqns. (66) to (71 ) become uncoupled and 
can be solved indepentJent of each other. For beams of mono- 
axis symme'cric cross-soctions, iqn. (133) holds and the 
Sqn, (66)to-(67) become uncoupled and the Fqn, (68) to (71) 
remain coupled. 

If the beam cross section hns no axis of symmetry 
, tJhian ai i tka Acinsti nns (661to(71) remain coupled and have to 
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^B-y-2 


0.04 0,06 0.08 


(a) Pinned -Pinned Beam 
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(b) Fixed -Fixed Beam 

% Decree in Frecpjencies of on. yfiN^qual Channel 
Secikm Mmm ap FtvK^ons Plj&fev. 









3,3,1 CuntilevGr oC box cross-r.ection 


ilv*re a .led cantilevjr ol I'o.' cross oectioii, 

which is doubly symmetric, Fi^.O, used by [7 ] is 

consiuered , Ihe dimensions md tho section iDropertics of the 
beam are given in Fig. 9 , derivation of these section properties 
is given in Appendix III. as this cross section is doubly 
symmetric, governing Eqns. ( 66 ) to (71) are uncoupled. The 
static analysis of beams in bending is v/ell known in the 
literature. Hence the static analysis of thin-walled beams 
in torsion, using warping -torsion theory is given in this work. 
Thus governing equations 'ire (70)to(71). It is noted that 

Che formulation of Chapter 2 does not include the strain 
energy due to St*Venant shear stresses, Chen and Hu [ 7 ] » 

v/hich is 

1 

.... (140) 


V = GJ^ f dz 
o 


where 


\ dA 


.... (1 41 ) 


is the 3t Venant torsional constant corresponding to 
St Venant shear stress. 

Inclusion of this term will give rise to the 
following term in ^qn. (70) 


sv 


2 Is i 


.... (14Z) 
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Fig. 9 Thin -Walled Cantilever of Box Cross-Section. 
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J . 


o 


.... (143) 


Thus the results in this section are based on i^qns . (70)t»(71) 
with Hqn. ( 142 ) included in Son. (70). Loading is uniformly 
distributed torque, m.j = 0.0216. The bound .ry canditions .ire 


f 


Fixed end 


Free end 


= 0 , f = 0 



§ § 

SeL Sv^ 


= 0 


(144) 


.Results for angle of tv\'ist (Q^)* v/arping measure ( 0 ) 
and bimoment (M g) are given in Fig. 10. Chen and Hu [ 7] 
results are also given and present insults match well. 



3.3.2 Tapered cantilever of I section : 


In this case a thin-wallod t''*pered c'-ntilevcr of 
I section used by [S ] > chosen. This section is a 

doubly symmetric cross-section. Thus for torsion analysis 
governing equations are Sqns . (70)to(71) with Eqn. (132) and 
boundary conditions are Eqn. (I44j. The dimensions and section 
properties are also given in Fig. '11. Squotions (70)to(71), 
which can be used for variable cross-sections, are solved for 
a end torque of = 224..65E03 (1 kN.m) . Ten elements are 

taken in this case. 











20.0 


mr = 0.0216 



-2.801 



Fig.K) Variation of the Nodal Variables of a Conti 
lever of Rectangular Cross Section Along 
the Length of the Seofft. 




i Section A -A 



Fig.11 Thin-Walled Tapered Cantilever of 1 Section . 



87 


The re;.ul‘.-r> for of tv/ist (0,) and bimoment 

(Mq ) nre given in Fig. 12 as f .notion of taper (a= * 

For uniform cantilever (a=1) close form results are available 
in [2] present results match v.'cll with them . ‘"Results 

for angle of tv/ist obtained by j* 3 ] are -'Iso included, and 
present results match well with them. It may be noted that 
[ 8 ] divided the beam into elements ^long its length and 
then, a shell mid surface of the element is approximated by 
arbitrary triangular subelements. It is seen that bimoment is 
maximum at the fixed and it increases \/ith r , as expected. 


tvarping stresses (normal and shear) at the fixed and 
free ends are given in Table 8. 'f'hese are calculated from 
bimoments and tv/isting moments (at tho nodes) obtained from 
the elc-mont equ-tions and further using tho following 
relations [ 2 ] . 


0 




My(z) W(S) 






t I 




« • • « 


(145) 


where (S^) = . / u dA and Sg is the distance from the edge. 

For uniform cantilever ^a=1) close-form solution of 
[2] is also given in the table, present results agree well 



with them 



0zx1O 




Fig.l2 Angle of Twist and Bimoment Diagrams of 
a Tapered Cantilever of I Section . 


— OCXD O 

« * * r . 




Table 8 : Warping stresses in a Tapered Cantilever of I-seotion 
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3.3,3 iuono-axis s/mmetrlc be&m : 

For this study a thin-walled ^ 

j uxLi waxiea toeam of semicircular 

section, fia. 5, Is chosen. For this mono-cxlsymmetric beam, 

as said earlier the govornlni equations -.ro liqns. (68)to{71) 

'.vlth (133). The hound-ry conditions -.re Eejn. (135). 

Biis beam is subjected to a unilcrmly distributed load., 

qy =■ -1.0. It Is solved lor pinned- pinned, .Uxed-fixed and 

cantilever coses, unc.^lts for displacement in the y^irection 

(U„). ansle of twist (e^) and warping n,easure(0 )are given in 
B'igs, 13 to 15. 


rigure 13 gives the results for pinned -pinned beam. 
Displacement and angle of twist are maximum, and warping measur 
is zero at the center of beam as expected. Maximum displacement 
io 0,01 uu1 3 cornfiared with elementary strength of material 
value of 0.013021, 

‘-Results of "^ig. 14 and 15 for fixed -fixed and 
cantilever beams are also as expected. Maximum displacements 
are 0,0028963 and 0.125168 comparcii to elementary strength of 
material values of o. 0026041 and 0.125. 

3.3.4 Unsymmetric section beam : 

For this, a thin-v/alled beam of unsyinmetric channel 
cross section. Fig. 7, is chosen. As the cross section is 
unsymmetric, all the equations (66)tD(71) are coupled. These 
coupled equations are solved for pinned -pinned, fixed-fixed, 
and cantilever cases. The bourdary conditions are same as in 




Fig.13 Variation of the Nodd Variables of a Pinned -Pinned 
Beam of Semicircular Section Along the Length of the 
Beam. 






Fig.14 Variation of the Nodal Variables of a Fixed -Fixed 
Beam of Semicircular Section Along the Length of 

Beam. 
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Fig. 15 Variation of the Nodal Variables of a Cantilever 
of Semicircular Section Along the Length the Beam. 
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Eon, (138). These bourns r;ro subjected to u uniformly 
distributed load, = -1.0 


The results for displacements (U^, U ), angle of 

y 

twist (G^) and warping measure ( f ) are given in Figs. 16 
to 18 xor these beams. Results are as expected. Displacements 
less than the displacements U,^. llaximun displacements 


Uy are 0.013 j> 72, 0.C029552 and 0.12388 compared to the 
elementary strength of raatoriol valuas of 0.013021, 0.0026041 


and C,123 for iiinned-piunud, fix'-;ci 


Ixo'-l -ind contilPV<'‘r cases 


respectively , 


3.A 


ITUCIXI 


APALYblO 


4 


For a column with doubly symmetric cross section, the 
shear center B and the centroid 0 coincide. 

Thus 

Kg - Yg = 0 .... (146) 

Hence the equations (93)tD(95) become uncoupled and 
can be solved independent of each other. For a column with 
raonoaxisyrnmetric cross section, shear center lies on the axis 
of symmetry. If the x-axis is the axis of symmetry, Fig, 5, 
then 

Yg = 0 

Hence the Er^n. (93) becomes uncoupled from the rest 
of the equations, ^qns. (94)b(95) remain coupled. This is 
known as double coupling of flexural buckling about y-axis 
and the torsional buckling. ■ 





Fig.16 Variation of the Nodal Variables of a Phned-Pinned 
Beam of Unequal Channel Section Along the Length 
of the Beam. 
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In the CHS 
shear center docs 


e of ci column -yith no c^xis of symmetry, 
not lie on any axis an.il hence the governinj, 


equations, inns. (95)t»(95), 
triple coupling of flexural 
^nci torsionai buci:ling a’oou 


remain coupled. This is knovm as 
buckliiigs about x and y-axes 
t 2 -axis. Bouiur ry conditions are 


Mins. '.97 ,^(99). 


i^or buckling analysis, unsynmetric channel section, 
Fig., 7, is chosen. Crdtic"! lo'^d-s are obtained for pinned- 
pinnod, fixed-fixed and c.'.ntilover columns as function of 
^®xx given in ^igs. 19 to 21 respectively. Close form 

results for pinned -pinned and fixed -fixed beams of [l?] 
are also given in the figures. Present results match well 
v/ith them. 


The critical buckling load F increases with increase ol 
as expected. It is observed that this increase is not 
uniform. In some modes the increase is Shnli and in some modes 
it id vary large. This is because in couploti buckling, there 
will bo some modes v/hich are flexur'.a dominant and some other 
modes are torsion dominant, dtudy ol mode shapes, for pinned - 
pinned column shows that the first' and third modes are torsion 
dominant, second and fourth are flexure about y-axis dominant, 
fifth mode is flexure about x-axis dominant. 



0.03 0.04 0.05 0.06 0.07 0.08 

Rexx ► 

Fig. 19 Critical Buckling Loads of a Pinned-Pinned 
Column of Unequal Channel Section. 
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p »| — P X 



0,03 0.04 0.05 0.06 0.07 0.08 


Rexx ► 

Fig. 20 Critical Buckling Loads of a Fixed-Fixed 
Column of Unequal Channel Section. 




PxlO 
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Rexx ► 

Fig. 21 Critical Buckling Loads of a Cantilever Column 
of Unequal Channel Section. 
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i.'ie -li: rnoncioi'!;;] p. 'raniu torr'. .for thii; socti 


on 


are 


= C. 07-37 ; ROyy ^ 0.0282 ; = 0.001782 


Parameter Re is £:renter than Re„„, and in turn Re is 

JUt YY’ yy 

greater thaji There Tore increase in critical loads with 

increase in is l-srciest in flexure deminant modes about 

::-.a:ris and smallest in torsion dominant modes. Similar trends 
are observed in the case of fixed-fixed and cantilever 
columns . 
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3.5 PLANE PRAMS : 

In this section results for the static and free 
vibration analyses of two plane frames, a balcony frame 
and a chassis frame subjected to transverse loading. Fig, 22 
and 24 respectively, are given. These frames lie in the xz 
plane. Beams of channel cross-section, shown in Fig, 23, 
are used in both these frames . The dimension and its section 
properties are also given in this figure. 

3.5.1 Static analysis : 

3. 5. 1.1 Balcony frame ; 

In this case a balcony frame, Fig, 22, is considered. 
This frame is formed by joining three thin-walled beams of 
channel cross-section. Fig, 23. A concentrated load of IkN 
is acting at the mid point of the longer beam. It is acting at 
the shear centre , The frame is divided into eight finite 
elements . Firstly the element stiffness and load matrices are 
calculated in local coordinates and then transformed to the 
global coordinates using the transformation matrix. These are 
assembled for the frame , After applying the boundary conditions 
the displacements are obtained* Stresses are calculated from 
Eqns. (123) to (124), after obtaining the elemental forces 
and moments by back substitution from Eqn. (122). 
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Fig.22 Balcony Frame. 


y 



^ $ A ** S 

Ixx s 2.4943 X 10 m, lyys 0 4787 xlO , A s 1,66 x 10 x 

Xb s -.33mm Vfes 0.0 mm Id s 42.0 xIO’V, Us 1.053 xio’m® 
lo* 2.973 X 10 m*, Ib= 4.781x10 m^ E = 200 GPa 
Gs75GPa, f = 78B0Kg/m^ a* = 20.506 mm 


Fig, 23 Thin-4/Valled Channel Cross-Section. 
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^■19.24 Chassis Frame. 
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Displacement, are given m Tabl» 9 a^d .t 
Table 10. Firet result, are obtained b 

of warping, i.e. Euler be ^ "«els=ting the effect 

b. i.e. muler beam elements are used Tb= b ■ 

equations for this case are f . * ^ 

qns, (104) and (105) with thp 

omission of the term with warnln. • 

n warping rigidity (El ) in 

Eqn. (105). In this case linear intern 1 . ““ 

used fnr +K^ 4- • interpolation function is 

used for the twisting anele o m 4. 

obtained with 

including the effect of warping i a tvi- 

P g» i.e. thin-walled beam 

elements are used. The results of T 1 r 1 v. . 

s+iff ^ obtained by direct 

stiffness method are also given in the Table 9. 

The present results for displacements and stresses 
based on Euler beam match well „ith Ma] . ^e displacements 
b»3ed on thin walled beam elements (i.e. inoluding warping) are 
quite Close to those based on Euler beam elements, but stresses 
dlTfer quite a bit. ^ere is not much change ip displacements 
as the load acts through shear centre and causes no torsion. 

On the top flange, sectorial coordinate u is positive 
maximum at A, zero at C and negative maximum at D , Fig. 23 
Shear stresses are maximum at the C point. The stresses are 
calculated at these three points and given in Table 10. Thus 
warping effect changes the stresses quite a bit. 
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3,5,1 .2 Chitssls frame ; 


For the analysis of a chassis frame, Fig. 24, is 
considered. The dimensions of the frame are given Fig. 24. 

The frame is made up of seven thin-walled beams of channel 
section Fig. 23, The dimensions and section properties of the 
channel section are given in this figure. The frame is 
divided into 36 elements with 33 nodes . The frame is simply 
supported at the nodes 2,8»1^ this present analysis 

two configurations of the chassis frame are used.Tv?© 
different types of loading are considered. Results are obtained 
without and with warping effects. The following notation is 
used to identify the different loadings and configurations. 

In all the cases the loads (along y-axis) act through the 
centroid of .the cross-section. 


B1 - Unlformily distributed loading of 8 kH/m on 

oireumferenoe of the frame only . All the circum- 
ferential Channels are facing inside. Beam 3-- 19 

parallel to beam 1 -Zl; I'bama 5 - 17 and 7 - 15 
are parallel to beam 9-13. Effect of v^rplng is 

neglected in this case. 

47 -heams is same as in B1 . Vfarping 
WB1 - Configuration of yearns 

effects are included* 


is 
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BG1 - Loading is same as in case. Longer beams face 
inside. Local x-axes of all shorter beam are in 
the direction of global negative z-axis. Warping is 
neglected . 


WBC1 - Configuration and loadings are same as in BC1 , 
but warping is included. 

Cases B2, WB2, BC2 and WBC2 have same configuration 
of beams and warping effects of B1 , WB1 , BC1 and WBC1 respectively , 
but uniform loading of 8 kN/m is on all the beams. 


Finite element equations for the whole frame are 
obtained as explained in the case of balcony frame ; displacements 
and stresses are obtHined. The nodal displacements (of shear 
center) and rotations for B1, WB1 , BC1 and WBC1 are given in 
Table 11. Studying the displacements for B1 and WB1 . 


It is observed that there is no change in the displace 
ments along the nodes of the longer beams 1 - 9 and 15 - 21, 
a bit of Change in outer shorter beam 9 " 13 and some change in the 
other shorter beams. Thus the change in the displacements, vfcen 
warping is included, is observed along the shorter beams . 

Comparing the displacements of B1 and BC1 oases there is a 
Significant change in displacements along all the beams. This is 

because in the BC1 case all shorter beams are in the same 

. i-nramred to B1 case where 

direction and given rlie to more twist compar 

twisting gets cancelled because of four shorter beams face 


each other. 



Table 11 : Nodal displacements and rotations, of a chassis frame 
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Con,p..rlnr, the displacement 

01 BC 1 and WRn -*4- » 

observed that there is no change 

longer beams and the change in th . along the 

beams S - 13,7-15.5-17 and th® shorter 

as that of B1 and WB 1 oases. The d- 

1 91 4 « ^^Plscements of shorter beam 

1-^21 increase in contrast to thein ^ 

decrease between 

B1 and WB1 cases. 

The total normal stresses 4 .^. 

the points A and D, and 
shear stresses at point C are given ^ -n 

in Table 12 for B1 , WB 1 
BC1 and WBC1 cases. The stresses .. 

"■or the cases B2, WB2, BC2 

and WBC2 are given in Table 13, 

Referring to the Tables 12 .^4 

od I 9 , it is observed that 
the stresses change when the effect 4 . 

of warping is considered. 

This change is very significant and c«n 4 

can not be Ignored. 

Also the stresses in BC and WRr- 

WBC cases are in general 

higher than those in cases of B and vn . 

WB, as in BC and WBC cases 

twisting moments of shorter beams do nm 

0 noL oppose each other, 

3.5,2 Free Vibration . 

For the free vlhretlon analy,i 3 of the ohaeels frame, 
the Chassis frame „f Fig. 2 , 1 , 

supported at the nodes 2,8.14 and 2o. elemental mass end 
stiffness matrices are obtained for each element in the local 
coordinates using Eqn. ( 127 ). Then these elemental equations sre 



Table 12 : Stresses in the chassis frame for loading case 
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transformed to the global coordinates, using the transformation 
matrix of Eqn. (115). nssembling for all the elements in the 
global axes, using Eqn. (13o) and applying the boundary 
conditions for simply supports , Eqn, (118) , the Eqn. (131) 
is solved for the frequencies of vibration. The first six 
frequencies of vibration are given in Table 1^. 

First the frequencies are obtained, neglecting the 
effects of rotary inertia and warping. The basic equations for 
this case are Eqn, (125) and (126) with the ommission of rotary 
inertia and warping terms. In this case linear polynomial 
approximation is made for 9 . Next results are obtained including 
the effects of rotary inertia and warping. In this case the 
basic equations are Eqns, (125) and (126), Refering to the 
Table 1 4 it is observed that the frequencies of vibrations 
decrease with the inclusion of rotary inertia and warping, as 
expected. f>ecrease in the fifth frequency is very small. It 
seems fifth mode is twisting dominant mode. 
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Table 14 : ‘^'requencies of vibration of a chassis frame. 


(rad /sec) 


Mode 

No. 

Excluding warping and 
rotary inertia 

Including warping and 
rotary inertia 

1 

289.80 

279.65 

2 

386.31 

360.80 

3 

412.62 

364.32 

4 

418.57 

369.39 

5 

423.24 

416.10 


6 


606,69 


524.86 
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CHAPTER- IV 


CONCLUSIONS 


The following conclusions are drawn from the 

present work, 

1 , ■The application of finite element method to the 

coupled equations of free vibration, static and 
buckling of beams and frames is very straight forward, 
even when the effects of warping, rotary inertia and 
shear deformation are taken into account. To find the 
exact solution in these cases is very difficult. The 
exact solution for frames with various types of beams 
does not seem to be possible. 

2. When the effect of shear deformation is included, 

accurate results are obtained using six coupled second 
order differential equations with cubic polynomial 
approximation as geometrical and natural boundary 
conditions are satisfied exactly. 

5. In the case of coupled vlhrations, v/hen the effects of 

rotary inertia and shear deformation are considered, 
the frequencies deorease. The percentage decrease in 
frequencies of higher modes is not always more «ian the 
percentage decrease in the lower modes, in contrast to 
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the uncoupled vibrations. This is because in coupled 
vibrations, there is an interaction between the bending 
inodes and torsional inodes hence some modes are bending 
dominant modes and some are torsion dominant modes. 

Since the effects of rotary inertia and shear are more 
for bending than for torsion, the percentage decrease 
in frequencies of bending dominant modes is more and that 
of torsion dominant modes is less. 

4, Present formulation can take care of tapered beams in 
a straight forward way, 

5 . The normal and shear stresses change quite significantly 
when the warping effects are considered. The warping 
stresses can be as large as or even larger than the 
bending stresses, hence can not be ignored. In frames 
stresses will change with change of configuration of 


the beams. 
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APPEIDIX -I 


DI3PI/tCEMSNTS OF A POINT OK THE CROSS-SECTION 


11: is assumed that the cross-section of the thin- 
walled beam retaias its shape when it bends and twists [ 2 ]. 
Thus the cross-section moves as a rigid body and such a 
movement can be considered as a rot'.tion about a point 
(Vlasovas theorem). 

Now a geometrical relationship which exists between 
displacements when a cross-section is rotated about an 
arbitrary point * Fig. A.1, is established in this 
appendix. 


The in -plane displacements : 


The displacement u (z) and u (z) of a point A on the 

cross-section due to a rigid -body rotation are derived in 

tei'ms of the displacements u (z) and u (z) of the origin 0 

° - ^o 

(Fig. A.1).' Points 0 and A move to O' and A respectively. 
Body rotates about point C through angle ©^. 

It is seen from the Fig. A.1 that 

FG = OE + (/H-g"h'- of 


or 




A 


... (A-1) 


and 


or 


Af'*' = FF'+ Hh'+ G a' -FA 




.... (A— 2 ) 
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nSLiumin,; anf.le of x^utation ©„ as small 
sin and cos 1), one gets 


(i.e, first order theory, 


= u 

X X 

- y.' 

Z 

.... (A-3) 

M 0 



i a U 

yo 

CD 

.... (A-4) 


Similarly displacement of any other point B can be expressed 
as 

%=%-yB®z .... (A-5) 

u. ,=» u„ + Xo 0_ .... (A-6) 

Yb Yq ^ z 

Thus displacements of point A on the body can be 
expressed in terras of the displacement of any other point B, 
as 


in. 

.... (A-7) 

“y.= ‘Vb'" «Z 

(A-8) 


The above displacements ; u etc. are in the 
directions of the coordinate axes. It is desirable to know 
the displacements in the tangential and normal directions of 
the section. A section VA is shown in Fig. A-2. Tanget at A is 
making an angle a with the x-axes . .Displacement in the 
tangential direction is v and in the normal direction is w. 

It may be noted that both these are in-plane displacements. 

It is seen from the Fig, krZ that 
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v(z,s) = cos a+u sin a 
A J A 

v/(2,3) = c +u cos a **** 

‘A 


U.'ji nc 

V(2,S) =U 
W(2,S) = 


Mjni 




cos 


(A-7) and (A-8), these equations become 
a+ Uy^sin a + ((x,,~Xg)sin a-(y^-yg)cos a)©^ 


•u^sin c.+Uy^cos a+((x^-Xg)cos a+(y^-y^)sin a)©, 


A 


.... (A-10) ■ 

From Fi", a-2, it is seen that perpendicular distances 

(Pjj and q^) from B to the tangent and normal can be written as 



BB' = EA^ - B a' 



or 

Pb «(x^-Xg)sin a- 

(yA-yn) a 

.... (A-11) 

Similarly BA" = BE+b''a 



or 


(yA-yB^sin a 

.... (A-12) 

Also 

dx = ds cos a 



or 

cos c=® - X 


(A-13) 

Similarly sin a = ^ = y 


# • • « (A— 1 4) 


Here 'bar (-) denotes differentiation with respect to 
s .coordinate . 

Using these Eqns. (A-11) to (A-14), Eqns, (A-10) 


become 

v(Z,s) = Pb ®Z 

' -"xg *z 


(A-15) 


(A-16) 


# • 
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The out~o r>T)l,-ine di., placement for 

iii — Oi)en proLile sections : 

_'he shear strain 7 is 

2S 

'zs ~ *’ 

■■Vhere Is the shear strain due to warpinj and r is the 

s 

shear strain due to daint Venant torsion. For open profiles 7 ^ 
is of seconaary nature and is assumed to be zero at the 
middle surface of the profile (Vlasov’s theorem). The 
Saint Venant shear strain 7^ is also zero according to Saint 
Venant torsion theory for the open profiles, thus 

2 S 3s 3z ^ *••• (A-17) 

Here u (without any subscript) is the out of plane displacement 
i.e, longitudinal displacement. 

By integrating the Eqn.(A-17) between an arbitrary starting 
point V and point A on the profile distance s from V, Fig. a- 2 
and .substituting v from Eqn. (a- 15), the following expression 
for the displacement in the longitudinal, direction is 
obtained . 

u(z,s)= Uy.(2) ~ x(0) )-Uy^(z) (y(s)-y(O) ) 

-©^(z) ...• (A-18) 

where (<Jg(s) is the warping function and is 

s 

Wg(s) = / Pg(s) ds .... (A-.19) 

ri 
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Ubv..ju.,]y iu-o x(o) mid y(o) are cunstani, Cor a given 
V.-Iuc o. ... ..e can, therefore, Incorporate their products 
v.'itn conctant u,,(z) into one constant [5 ] and rev/rite the 
above euuution (A-le) into a general form : 

UU,S)= U^3(2)-.u;^(z) x(s)-u;^(z) y(s)- Q^Cz) c^(s) 


(A- 20 ) 

-displacements Uj^jU^u^and rotation 0 ^ can be considered 
to have two parts . One part with subscript b denotes the 
part of solution in which shear strain has been neglected and 

the subscript s denotes the contribution of the shear strain. 
Thus 


% b' % s 

""ypr % b" s 


0 = 0 
b ^ 


= U + u 
® % ®s 


(A-21) 


Using the Sqn. (A-21), hqn. (A- 20 ) becomes 


U(2,S)= U^(z.s)tUg(3,s) = Cu3^-ui^ > 

D D 

.... (A-22) 
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Where u^(z,c) is the out- of -plane displacement at point A 
neglec tir.js, the s.^ear strain and Ug(z,s) is the out-of-plane 
displacement at point A due to shear strain. Displacement 

sraoll compared with u^(2,s) and can be neglected 
16 J , Thus 

u(2,s)= u^(e,s)= Ug^ -u' x(s) V y(s)-. ©' w^(s) 

h c b b D ^ 


(A -23) 


u!. and u/. 


"’'V“ 


3u. 


b 

b 


in £qn, (A-23) are 


8u. 


u 


b 


d2 


and u 


’^B b 


'^B b 


3z 


From Fig. (A-3) it is seen that above derivatives can 
be written as 


u:. « 0.. and u.^. =-6. 


^ b 


^B b 


.... (A-24) 


(z) in Eqn. (A-23) is measure of warping and now onwards 
^b 

it will be denoted as 0(z). 

Now onwards u_ in •t^qn, (A-22) -is denoted as u for 
^b 

the sake of simplicity. Thus Eqn. (A-23) becomes 

u(z,s) = u^{z) - ©y x(s)+ 0^ y( 5 )- 0 (z)(^(s) 


• » • • 


(A-25) 
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y 



Fig.A-3 Variations of the Displ<xement Ux and Uy along 
the Beam Axis. 
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inhere ’i’j_(z) is the torsional function, independent of the 
siciteri^X and external loading of the heam. 

Hence 

_ au(z,s) 3v(z,s) f^( 2 ) 

'^zs~ Ts ~a 2 “ “T'Cs) •*“ (A-30) 

Now integrating Eqn, (A-30) between an arbitrary 
starting point V and point A on the profile, distance s from V, 
and substituting v from the Eqn. (A-15) and using the Eqn. (A-21 ) 
and neglecting the out of plane displacement at A due to 
shear and using Eqn. (A-24) the following expression for the 
displacement in the longitudinal i.e, z direction is obtained, 



.... (A-31) 

s f.(z) 

Here / (p 2 (s)~ is defined as the warping function 

to>g(s) for the closed profile 

s ^^(z) 

o)g(s) *» / (Pg{s) - ),, ds .... (A-32) 

Now Eqn, (A-31) becomes, 
u(z,s)» Ug^(z)- ©y(z) x(s)+ ©^(z)y(s)- 0(z)wg(s) 

.... (A-33) 

summarizing, out of plane displacement u and in -plane 
tangential and normal displacements v and w for thin walled 


section are 
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where 


Here 

point 


uC2,s) = Ug- 9yX+ e^y . 
v(z,s) “ X +Uy y + p 

v/(2,s) = -u^ y+ Uy X + q 

*!• 

C* 

“ / p ds for open profiles 

o 

s f . 

w = ; (p- :^)ds for close profiles 


.... (A-34) 


(A-35) 


^X* ®X» ®y» ®2» P» ^ W all are of 

chosen. 
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•■■r?Sl€)lA^ 7T 

9MIV^-TI0* ^0F PniK'CIPAL POLS AI'E) PRINCIPAL ORIGINS 

Consider v. thin-vvallod cross -sect! on , Fig. B-1 [3 ] 
^fhich io hriving tv/o contour coordinate systems s and s^ 
and coi*ru3ponding origins 0 and 0*. 

p(^^»yp) is tne pole oT the cross section in the contour 
coordinate- system s and 

P is the pole of the cross section in the contour 

coordinate system 

The perpendicular distance from the point p to the tangent 
at any point A(x,y) on the cross section of the heam is r and 
tho perpendicular distance to the normal at k is q • 

Similarly the perpendicular distance from the point p to the 
tangent at any point A(x,y) on the cross section of the "beam 
is r and the 'perpendicular distance to the normal at A is q . 

The coordinate of 0 in the s contour coordinate 
system is s^*. 

There will be two contour v/arping functions («>(s) and 
to (s ) associated with the points pAO and p AO respectively . 

By the definition of contour v/arping function 

to(0) = 0 

to*(0) = 0 


.... (B-1) 
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Since tho 
position 
r chances 


perpendicular Oir.truices r 
of the pole, if the pole i 
to r^. 


and n are dependent on the 
s changed from p to p*, 


It is seen from the Fig. b- 1 that 
= r(s) - Apx- 

“ r(s) ~ BF*+BA = r(s)-BP^+CD 
r*(s*) = r(s) -(Xp# -Xp)sinG +(y^^.-yp)cos Q 
n*(s) = q(s) - .'.P 

= q(s) - AD - DP = q(s) - BC - Op 

0 0*0 \D J 

q*(gf)a q(s)- (Xp*-Xp)cos 0 _(y^^-y^^)sin © 
s* 

since w* a / r* ds**- 
o 

a change in w* is .... (B -4) 

d(i)* a r* ds* 


Since a change in s* , is equal to change in s i.e, 
ds* = ds , Fig, B-1 , 

dw* = r* ds .... (B-5) 

substituting for r* from Sqn. (B^2) in (B-5) one gets, 
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du'' “ r(G) do -(Xj^^-x^)sin Q ds:-(y ,(.-y ) qqs © ds 

.... (B-6) 

..e loiov^ th'it sin © ds = dy , cos © ds = dx [ 2 ] 

also r(s) ds = dcu(s) .... (b- 7) 

Substituting Eqn. (B-7) in Eqn. (B-6) one gets 

« - (Xp*-Xp) dy + (yp*-yp) dx .... (B-8) 

Integrating Eqn. (B-8) w.r.t s between the limits s^* and s 
i.G, bct\'/een the origin 0* and point A one ge'ts 

«* = <^-%*--(Xp4t-Xp)(y-yQ^)+ (yp*-yp) (x-v^ 


whore 


.... (3-10) 


w 4. is the value of w at 0* and it gives the difference in 

o* 

the origins 0 and 0*^ (^"'iS* 3-1). 

Now the properties of the cross-section w.r.t. the contour 
coordinate system s* can be expressed as 



/ t# dA 

A 


Vx” 

^ 0*y dA 

A 

.... (B-11) 

I = 

w*x dA 


t«3*y 

A 
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Un .“.iil-.'.t:! f.' trio v .lue 
nnd siinrlyfyiiv-, one ^'ets, 


of from Sqn. (:,._9) (B-11) 






( f ’=‘ 1 “ -^o* i dA) 


!„*= A -(Xp, -’'p ) (ix-yo» '• ) * (yp» -ypX ly-X;,, A ) 


.... (B-12) 

^o)*x=‘ { <‘'*yiA " ^ (“y-“o*y-(V’'p^<y^- 5 'c>^y)"(V-yp)(xy-Xp*y))dA 


=■ MYdA-.p. /^ydA-(Xp*-Xp)(/^y 2 dA-yp» )*(yp,-yp) 


( ^xydA- X „ /ydA) 
n A 


-<.>*X =lBX-“<fIx-<'x,»-=fn) (Ixx-yo*^X^ + (yn»-yr.) 


'p' ^xy 


...■.(B-13) 

similarly 

(B-14) 
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The ^iv.= the to the sectional prooar tl -s 

v..ho-, rov a,.-: ■■■-.;i^'-reohon,:ed.IntheaboveequatioL.. 

x,y axes, are arbitrary. , 

'..hen the axes i,y are the principal axes .aod the origin is 
located .at the centroid of the cross-section, one has 

I..as / ydA = 0 


I = / xdA = 0 


(B-15) 


^XV~ ! ^ 

j A 


By raakin^ use of kqn. (B-15) the Jiens. (3-12)to (B-15) 
reduce to 



.... (B-16) 

■W*X * ^ox“ ^xx 

• . » . (B —1 7 ) 

6)*y ” ^oy'‘‘^^p*"’^p^^yy 

.... (B-18) 


Eons* (B-17) and (B-18) are independent of the position of 
the contour origin* Prom Eqns, (B-17) and (B-18) one can see 
that ^ 

^hen x_* ^ ^rv •••• Cb-19) 

p p ■^xx 

I S' 0 

, \ 
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y 



y. 


yy 




I 


w*y 


= 0 


The pole P*(Xp*,yp^) where the <iuantities and 

are becoming zero is called the principal pole 
(cenzcr o£ shear) . 

Hence when the pole of the cross-section is the principal 
pole (center of shear) then 




1 w - 0 


.... (B-21) 


From Scjn. (3-16) one knows that 

when the two poles p and p* are placed at the principal pole 


and 


(i) 

Wq* * r“ 


(B-22) 


then 

=- 0 


The position of the origin 0^ vdiich gives the value 
of Wq* (Eqn. B-22) is said to be at the principal position 
,i.e.!o» is the principal origin of the contour coordinate 

system* 
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iicncc '••lior. the origin (starting; point) of the contour is 
the principal contoux- origin then 

= 0 .... (B-23) 

The s* contour coordinate system is the irincipal contour 
coox'.’-iJ'i^'te system and all the section properties with a * symbol 
aro C'.lculated with reference to this syrrceni. 







APPENDIX- III 


C,1 Section Properties of a channel section : 

A channel section of breadth b and depth h, Fig. C-1 
is considered. The thickness of the web is and that of 
flange is t^.. In calculating the warping function U)(s) of the 
cross section, clock wise sense is taken as positive . 

Consider the pole and starting point to be at the 
point D, Fig. C.1(a). For finding the warping function of 
this open profile the following equation is used, 

w(s) =« J p ds (C-1) 

s 

where p is the perpendicular distance from the pole to the 
targent to the middle line of the profile at the point of 
consideration. 

At the starting point D, s = 0. Since p is zero for 
both DA and DE the warping function w(s) is zero from E to A. 
From E to a point lying at a distance s^ to the right of E is 

,,( 3 .) . - hds= -hs, (C-2) 

‘ 0 

Here the minus sign appears because from E to F the 
sense of u is anticlock wise,i.^ negative. Thus from 
q(s) is a -linear function and at F it is given as 
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Refering to Fig, C,1(b) , where the starting print is at D 
and the pole is at B. By following the above procedure for 
calculating o>(s), one gets 03(s) as a linear function from D 
to A and at A it is 

“ab “ ^ (C-4) 

From D to E again it is linear and at E its value is 




EB 


== a. 


(C-5) 


When the integration is continued up to F, the value of w(s) 
becomes 

“fB “ ^ (C-6) 

When the starting point is at G , pole is at B, 

Fig, C,1(c), u(s) = Oat G. Integrating from G to D gives a 
anticlockwise rotation hence w(s) at D is 


Integration from D to A involves a clockwise rotation 
and hence <*>(s) at A is 




AB 


=-a h 

X 7 . 


. hb 


(C-8) 


At F o>(s) is same as at A but with a negative sign. 

While moving from D to A (i)(s) changes its sign from 
-ve to +ve through a point where m(s) = 0. This point at which 
03(s) = 0 is at a distance of a„ from the web. At this point 



* paq-BinoXBo sq ubo saTqjadojcd uoTqoas 
iqo sqq AxaBipniTs 'uoxqoes snoTASjd eqq ut ‘uoxqoas x^tiuBqo aoj 
'qxaosap si3 paqBxnoXBO aq ubo (s)co uoxqounj SuxdjBM aqq */. ‘Stj 
‘ uopqoas x^uuBqo x^Bnbaun jo msaq paxx^.'^-'^mq -JO^ 

'"uoxqoas XB'i^tsiun ub jo saxqaadojd uoxqoag *Z -0 


q^q qg+Jq qZ+^q^^ 

(ZU-0) ^ — X — 


^q qz+^q q 

-,q 


®x 

^XJBXTtQTS 


(U- 0 ) 


^q q+^q qq ^ 

SB uaAxS sj ’^B aouBqsTP sqx 


(Ot- 0 ) 




AT '•I 

q ^z^■'^^ ^q^q -^q 


V 


3 uaAXS aq ubo puB pauxBqqo aq ubo yP (s)_c«) = AxJBXTinxs 


(6-0) 


Jo TI • lZ_ 


^q q 


.(Vq) 


. sp (s)n q ; = q -sp (s)n / = VP (s)ro / = (’‘p)'’i 

B-q 


<n. 


sx qaM aqq moaj b 
quBqsTp quqod sqqq qB ''i jo anpeA mnraxxBm aqq aouaH 

*^I quamoin XBX'^ioqoas 


XBOxqBqs jo anxBA ranmxxBm aABq xITa qi ‘0 = (s)i?) aaaqA 
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For this section the contour coordinate of the 
profile of the cross-section, s is equal to -x, -y and x for 
the top flange, web and the bottom flange respectively. 

Hence the section properties 1^^, 1^^, 1^^, Ip^ 

and of the unequal channel section are 

^Qo = dA = / (||)^dA = 9.03224 x 10"^ 

A A ® 

I^p = 5Ey dA = 0 

= ; p X dA = 0.283004 x 10”^ m^ (C-13) 
qr ^ 

I = J y^ dA = ^ dA = 5.80644 x 10“^ 

PP A A 

I = j p y dA = 0.36768 x 10“^ 

pr A 

and I = -■ p^dA = 0. 465625 x 10“^ 

A 

The non dimensionalised parameters are obtained 
using the properties given in Fig, 7, and Eqn, (37). 

C.3 Section Properties of a box sec t ion : 

Considering a thin-walled box section, of length 2a, 
and breadth 2b, Fig. 9 the torsional function T is 2 
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f = 



2x4abxt 

'4'(a+b) ’ 


2ab t 
Xs+'b ) 



2 ab 

Ta^FJ 


I 


rr 


/ dA = 2b 2at +2a^x 2bt 
A 


(C-14) 


I 


d 


. / dA = 
A 



l6aS^t 

(a+b) 


Because of the symmetry the starting points V lie at the mid 
point of each side. The value of w is same for all the corners, 
but they differ in signs. The value of w at a corner of the 
cross section can be given as 




. 2abN ,, 


a+b 

(C-15) 


Using this warping function, the sectorial moment of 
inertia, can ba obtained as. 


I 


(iXi) 


at 


j «^dA 
A 



(a+b) 


(C-16) 
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